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THE COMPATIBILITY OF THE SURVIVAL PLATEAUX 
HYPOTHESIS WITH GAUSSIAN DISTRIBUTION 
OF POPULATION 


B. DAVISON 


COMPUTATION CENTRE, PHYSICS DEPARTMENT 
UNIVERSITY OF TORONTO 


It is often assumed that a highly peaked Gaussian distribution of popu- 
lation characteristics would imply a similar distribution of the survival 
probability. Such interpretation is, however, rather arbitrary and unduly 
restrictive. In this article it is shown, on the contrary, that, where the 
characteristics depending on many genes are concerned, even a broad 
survival plateaux will result, under certain additional conditions, in a 
highly peaked nearly-Gaussian distribution of the population character- 
istics. 


In applying statistics to biological problems we often encounter 
the Gaussian Law, and we are naturally interested in the interpre- 
tation we can give to the parameters entering this law, in particular 
to the position of the peak. It is usually assumed that this peak 
would occur at a point which in itself possesses some special 
properties, representing, say, the point of greatest survival proba- 
bility, or such like. Such interpretation is suggested by the analogy 
with the way in which the Gaussian Law arises in physics (for 
instance, in the case of an instantaneous point source in heat 
transfer problems the heat distribution at any subsequent time is 
Gaussian with its peak at the particular point where the heat was 
originally applied). However, in the discussion of the distribution 
of properties of populations, the above interpretation of the posi- 
tion of the peak is by no means the only possible one. On the 
contrary, it can easily happen that the peak of the Gaussian dis- 
tribution will occur at a point which has no simple explanation. 
For instance, if, as is usual in biology, a large number of inde- 
pendent characteristics are involved, then even a wide ‘‘survival 
plateau’’ for the individual may result in a highly peaked Gaussian 
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distribution for the population. The position of the peak will then 
correspond simply to the center of the plateau. This fact has some 
very important implications, as will be indicated below. 

We consider the following simple example. Suppose we deal 
with a certain species and let z,, 7,,---, #, be certain character- 
istics in which one individual of the species may be different from 
another. For instance x, may be the wing-span of a butterfly, 
w,—the length of its antennae, x,—the number of green spots on 
the wings, and so on. We shall assume that the characteristics 
@,++.@, are chosen so that they are all mutually independent 
(uncorrelated), and each of them is measured in the appropriately 
chosen units; in particular the unit in which any 2, is measured 
stays in the appropriate relationship to the unit in which z,, say, 
is measured. 

Let p be the survival probability, for the individuals of the 
species in question, up to the age of producing offspring, and sup- 
pose that p has a constant value p = p, as long as z,, 2 
for the given individual satisfy the inequality 


gosto ty 


N 

>: (ep aj)? <R*, (1) 

j=l 
while if (1) is violated » =0. We shall assume, for the sake of 
simplicity, that the species under consideration is not sex dif- 
ferentiated, so that each individual has only one parent, and sup- 
pose that if the parent’s characteristics are hs Seas zy » the 
probability density for the offspring to have the particular set of 
characteristics 7;, 24, ..., &, is given by 


ftit=rly, (2) 
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where we heve put for brevity | r—r’ 


We inquire under these conditions into the distribution of the 
characteristics x, ... a, among the individuals of the species. If 
a stationary distribution is possible the density function P(@, «+. &y) 
will be given by the solution of the integral equation, 


Bey --.2y)=d fi. (a ..a%)f([r—1'|) da’ ...de,, (8) 
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where A is the survival probability for the individuals whose char- 
acteristics satisfy (1), times the mean number of offspring per 
parent, the latter being assumed independent of parent’s -char- 
acteristics. 

This is an eigenvalue problem, i.e., the solution is possible 
only for a certain discreet set of A-values. We may expect, how- 
ever, that the actual value of the survival probability in the sur- 
vival region (1) will depend upon the total population of this species. 
If this population becomes too numerous the multiplication of pred- 
ators and/or the limitations of the food supply will decrease the 
survival probability. The equilibrium conditions will be reached 
when A is equal to that eigenvalue of (3), for which the solution 
of the corresponding time dependent problem is stable; that is to 
say, for which any random deviation from the eigenfunction of (3) 
will decrease with time, rather than increase. Now the function 
FCT — r’|), though we do not know it in detail, should certainly 
be non-negative, so that (3) is an equation with a symmetric non- 
negative kernel. For the equations with such kernels it can be 
shown quite generally that only the solution corresponding to the 
lowest eigenvalue is stable in the sense explained above.* Thus 
the equilibrium conditions will be reached when 4 is equal to the 
lowest eivenvalue of (3). 

We now introduce a further simplifying assumption, namely, we 
assume that the mean-root-square change in characteristics from 
parent to offspring is very small compared with RF, i.e., 
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If this is satisfied and f(|r-—r’|) decreases sufficiently rapidly 
aaa Ce 
with increasing |r-—r’ yan eat oe 232 then it can be shown, 
much in the same way as it is done in the theory of Brownian mo- 
tion or in the neutron transport theory, that the solution of the 
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*This can also be seen, without resorting explicitly to the time de- 
pendent treatment, in the following way. The probability density pa, {ant @y) 
should certainly be non-negative. On the other hand it is known that in 
case of an equation with a symmetrical non-negative kernel only the 
solution corresponding to the lowest eigenvalue is non-negative, whereas 
those corresponding to any of the higher eigenvalues are bound to change 


the sign. 
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integral equation (3) can be approximated, except in the vicinity 
of the boundary of the survival region (i.e., except for re rice 
by the solution of the differential equation, 


N 2 
0“ : 
ae +K° d= 0 ’ (5) 
jel j 
under the boundary conditions 
$| = 0. (6) 
| r- r |= R 


0 


Here x? is some constant depending upon A and upon the form of 
the function f(|r-—r’|). The equation (5) with the boundary con- 
ditions (6) is again an eigenvalue problem and, if (4) holds, the 
lowest eigenvalue of (5) will be approximately equal to the value 
of x* corresponding to the lowest eigenvalue of (3). On the other 
hand the solution of (5) satisfying (6) and corresponding to the 
lowest eigenvalue can depend only upon = y but not upon the 
ratios (@, ~ Liq)/ bt Ci |. Hence equation (5) reduces to 
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We need the solution which is regular at |r- r, | =0. For any 
particular N such a solution is given by 


4(N-2) J 


ere) etate = const=(T =< Feh- 
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(x|r—r |), (8) 
where J() is the Bessel function. In biological applications, 
however, N will be very large indeed, so that what we are really 
interested in is the asymptotic behavior (8) for large N. This can 
be most easily obtained by disregarding (8), going back to (7), and 
solving it by a series expansion. Proceeding thus we can easily 
show that, for sufficiently large N, (a, ... z,) is approximately 
given by 


i Sgn 
b(w, -.. wy) = const + exp = >, @4-%;0)? : (9) 


But x”, as mentioned previously, should be the lowest eigenvalue 
of (7) under the boundary conditions (6). That is, in view of (8), 
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k is given by the lowest root of the equation 


Js cyagy (KR) = 0. (10) 


The roots of Bessel functions of an arbitrary order have been studied, 
and it is known that the lowest root of J, (y)=90 is given by 
y =n-(1+e(n)] where ¢«(n) > 0 (provided n is real and positive). 
Thus for large N we have 


N-2.f, (N-N) LN N 
ney eas ae 9 Op lte *) (12) 


and with this value of x the expression (9) becomes 


N w\]2 
d(a,-+.+@y) = const - exp ' aR? h (3) > @; = “| pees (12) 
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Now, (12) represents a highly peaked Gaussian, and its peak is 
situated at the point r=T, (i.e., 2, = yo) Cg = Pag +++ Ly = fv) 
whose only distinguishing characteristic is that it is the point 
furthest away from the boundary. 

Thus we see how a highly peaked distribution of the properties 
of a population, closely approximating the Gaussian Law, could 
arise without any corresponding peaks in the survival probability. 

This may throw some light on the controversy between what may 
be called the ‘‘survival plateau hypothesis’’ [according to which 
the survival law is roughly of the kind postulated in (1)] and what 
may be called the ‘‘survival peak hypothesis.’’ The general spec- 
ulations as to how evolution could have led to the present diversity 
of living beings, as well as observations of the actual conditions 
of survival and non-survival, tend to favor the ‘‘survival plateau 
hypothesis.’’ On the other hand the phrase ‘‘survival of the fittest”’ 
would seem to imply survival peaks rather than survival plateaux. 
Also the observed highly peaked Gaussian distribution of individuals 
within a species is often taken as a support for the survival peak 
hypothesis. The example we have considered shows, however, 
that the observed distributions are in fact incapable of lending any 
support to the survival peak hypothesis, and the case for the survi- 
val plateau hypothesis is correspondingly strengthened. 

The decision between the survival peak and the survival plateau 
hypothesis is also of some practical importance. Consider, for 
instance, the question whether helping a weakling to survive is 
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good or bad for the species. According to the survival peak hy- 
pothesis by helping the weaklings to survive we encourage devia- 
tions from the optimal set of characteristics and this in the long 
run should be detrimental to the species. On the other hand ac- 
cording to the survival plateau hypothesis by helping the weaklings 
to survive we increase Ff in (1) thus doing something to promote 
adaptability and to protect the species from overspecialization. 
That is, provided we deal with a survival plateau type of situation, 
helping the weaklings to survive should in the long run be bene- 


ficial to the species. 
RECEIVED 3-4-57 
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THE PHYSICAL RATIONALE OF BALLISTOCARDIOGRAPHY 


RUSSELL H. KESSELMAN* AND ARTHUR GRISHMAN 
CARDIOGRAPHIC LABORATORY 
THE MOUNT SINAI HOSPITAL, NEW YORK, N. Y. 


Ballistocardiography is a special case of the physical problem con- 
cerning the dynamics of interaction of two particles of varying mass. The 
internal forces result from cardiac activity and the external forces are 
those exerted by the bed supporting the body. Physical analysis makes 
it evident that the internal forces not only arise from moving blood but 
are a resultant of all motions taking place in body substance relative to 
a reference system which can be based ultimately on the skeleton. The 
fundamental problem in ballistocardiography is the determination of the 
magnitude of the internal forces as a function of time. 


Application of rigorous physical analysis is necessary for the 
further development of ballistocardiography. At present a variety 
of instruments are in use giving different types of records (Starr 
et al., 1939; Nickerson and Curtis, 1944; Dock et al., 1949). In- 
terpretations vary as to the significance of the deflections re- 
corded (Dock et al., 1953; Burger et al., 1953; Rappaport et Gis 
1953; von Wittern, 1953). Desirable is a common denominator in 
terms of the precise quantities characteristic of the dynamic physi- 
cal changes in the body produced by the physiological activity 
maintaining the circulation. It has not always been very clear just 
which physical property of the body’s motion is being measured in 
ballistocardiography. 

Instruments are available for conveniently measuring accelera- 
tion, velocities, and displacements alone or together (Pordy et al., 
1952; Smith and Brian, 1953). The character of these records de- 
pends both on what is taking place at the body generator and on 
the recording apparatus used. The stiffness coefficient, the damp- 


*National Heart Institute Trainee at time of this study. Present ad- 
dress: Woman’s Medical College of Pennsylvania, ‘Philadelphia, Pa. 
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ing coefficient, and the mass of the system, as well as the binding 
of the patient to the system, quite apart from anything inherent in 
the patient studied, affect the records obtained. For translational 
motion it is evident that the force exerted by the body on its sup- 
ports is the basic quantity. It does not matter what these supports 
may be; the force will be the same. That force is the essential 
quantity to be studied is a consequence of Newton’s second 
law of motion: change of motion is proportional to the force and 
takes place in the direction of the straight line in which the force 
acts. 

The clarification of these problems may be aided by a straight- 
forward analysis of the physics involved, although the physiologi- 
cal realities must be kept clearly in mind. Burger’s analysis (1953) 
is physically acceptable, but a number of assumptions question- 
able from a physiological standpoint are used. In our conception 
the blood (and some of the associated cardiovascular structures) 
on one hand and the remainder of the body on the other may be 
considered as the members of a system of two particles of variable 
mass. In the past some investigators (Dock et al., 1953) have 
considered that ballistocardiographic deflections resulted from 
specific stages in the movement of the blood itself in and out of 
the heart and into the great vessels. The mass, moving in relation 
to a reference system which can be represented by the skeletal 
structures, includes not only blood and cardiac muscle, but also 
adjoining organs set into motion by the contraction of the heart and 
the passage of the pulse wave. These include blood vessel walls, 
diaphragmatic muscle, pulmonary tissue, etc. However, the con- 
traction of the heart is the primary source of the forces producing 
this motion, wherever the site of the motion may be. 


Mathematical analysis. In order to simplify the discussion, the 
component of the translational motion of the body along a single 
coordinate axis (the y~axis) will be considered. Head-to-foot di- 
rections will be considered to occur along the y-axis. Let 2 be 
used as a subscript referring to moving blood and associated car- 
diovascular structures and let 1 be a subscript referring to the re- 
mainder of the body. Let m, be the mass of moving blood and other 
cardiovascular structures and let m, be the mass of the remainder 
of the body. Let y, be the position of the center of mass of m, and 
y, the position of the center of mass of m,. Let F,,, be the in- 
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ternal force acting on m,, F,,, the external force acting on m,, 


iy the internal force acting on m,, Hess the external force act- 


ing on m,, and ¢ time. 
By Newton’s second law of motion, 
d(m,dy,/dt)/dt=F,,,+ Fy .,5 (1) 
and 
d(m, dy, /dt)/dt=F,,,+F,,,- (2) 


Adding equations (1) and (2) we obtain 
d(m, dy,/dt)/dt + d(m, dy,/dt)/dt=F,., + F 


y 2ey? 
or 
d(m, dy,/dt + m, dy,/dt)/dt =F ,,, + Paty: (3) 
since 
ne ie F diy (4) 
by Newton’s third law of motion. 
Define 
M=m,+m™,, (5) 
and 
Y= (MY, + My Yq) (mM, + My)- (6) 


M is the total mass of the body and 7 is the position of the center 
of mass of the body as a whole. Then 

Md? 7/dt? = F (7) 
since F, ,, is zero, m, not being in contact with the exterior. Us- 
ing equations (2) and (4), we can then obtain the equation 


d(m, dy,/dt)/dt = - F,,- (8) 


ley? 


From equations (1) and (5) it follows that 
ad((M - m,) dy,/dt\/dt =F yi + Frey: (9) 
Substituting equation (8) into equation (9), 
d{(M ~ m,) dy,/dt\/dt = - d(m, dy,/dt)/dt + Fy .ys (10) 
or 


d(Mdy,/dt)/dt - d(m, dy,/at)/dt = - d(m, dy,/ dt)/ dt + Fyeyi Al) 
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Justification for neglect of d(m, dy,/dt)/dt may be arrived at by 
the following approach (Karreman, 1957). If it is assumed for this 
purpose that m, and m, do not vary with time, equation (1) yields: 


d?y,/dt? =(F,, +Fyo,)/m, 


Liy ley 
or 
d(m, dy,/dt)/dt =m, (Fy, + Fay) > 
*. d(m, dy,/dt)/dt =m, F,,,/m, + mF, .,/m, : 
Equation (8) gives: 
Fis, = ~ d(m, dy,/dt)/dt, 


which, when substituted into the above equation yields: 


oo ms 
d(m, dy,/dt)/dt = - — d(m, dy,/dt)/dt + — ee 
My mM, 
Substituting the derived expression into (11) results in: 


m 
d?y,/dt? - |- = d(m, dy,/dt)/dt +m, Fiey/m, be 
1 


— d(m, dy,/dt)/dt + F 


ley? 
which after rearrangement yields: 

d*y,/dt? = —(1+ m,/m,) d(m, dy,/dt)/dt + (1+ m/m,)F 4, + 
As m, << _m, this equation may be approximated by 


d?y,/dt? = ~ d(m, dy,/dt)/dt + F (12) 


ley* 


Equation (12) is the result of letting d(m, dy,/dt)/dt equal zero 
and of letting M be constant, the latter an obviously justifiable 
assumption. 


Using equation (8) and transposing in equation (12), 


Fj, = Md? y,/ dt? - (13) 


res. 
The force exerted by the dynamic activity of the cardiovascular 
system on the body is given by F’',,,. The determination of Fay as 
a function of time is a ‘paar goal of balistoen el ere cee 
The force F, , depends on the characteristics of the table and 
associated syaten) as regards the table’s mass, damping coefficient, 
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stiffness coefficient, as well as the character of the patient’s 
binding to the table used. In the direct body pick-up system the 
characteristics of the contact of the patient with the fixed table 
upon which he rests helps determine F’, ,,. Ideally the character- 
istics of F,,, should be such as to make F,,,, a function of time 
most readily available with a minimum of distortion. 

Consider the table types of ballistocardiograph and assume that 
the mass of the table is small and that the binding of the patient 
to the table is good (rigidly attached to table). Let & be the damp- 
ing coefficient of the patient-table system and & the stiffness co- 
efficient of the patient-table system in relation to its supports. 
Then we find 


F,.. = -Rdy,/dt - ky, , (14) 


ley 


or using equation (13) 


Fiiy = Md? y,/dt? + Rdy,/dt + ky,. (15) 


It is emphasized that FR and & are coefficients referring to m, and 
the relatively massless bed to which m, is rigidly bound. If the 
mass of the table is substantial, it is best to’ consider M as the 
sum of the weight of the body and the weight of the table. 

It follows that a determination by direct measurement of y,, 
dy,/dt, or fd y,/at? as functions of ¢, coupled with a knowledge of 
M, R, and k, would suffice to learn F’,,, as a function of time. 

The force F’,;, may be regarded as a complex wave form which 
can be expanded in terms of an infinite series of harmonic func- 
tions of the form 


F ity =), (Fyq 008 2nvt + Fy, sin 2nvt), (16) 


vel 


where a typical vth cosine term is F' cos 2nvi with v the frequency. 
Similarly 


Y, = a Yyo (17) 
vel 


in which y, is the displacement produced by the vth component 
term of F’,,,. A typical vth term of y, which corresponds to the 
vth term, F cos 2avt of F,,,, will be called y. The following dis- 
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cussion would apply in like manner if a typical sine term were 
used instead of a typical cosine term. Then 


F cos wt = Md? y/dt? + Rdy/dt + ky, (18) 
where wm =2mv. The steady-state solution of equation (18) is 
y = F cos (at - 4)/Vw2 R? + (k - Mw?)?, (19) 


where 
tan & = wR/(k - Mo”). (20) 


The second derivative of equation (19) is given by 
d*y/dt? = ~ F cos (wt — 4)/M VR2/w2 M? + (k/w2M -1)2. = (21) 


For ballistocardiography the choice of FR and & should be such 
as to make (/R2/w2 M2 + (k/w2M —1)2 as close to unity as pos- 
sible and to make wR/(k — Mw”) as close to zero as possible for 
the significant ranges of w. When these requirements are satisfied, 
measurement of d*y/dt? gives a result whose amplitude is inde- 
pendent of w. Expressing this mathematically 


D2 a2y,/dt? = - =)” F, cos a2, (22) 
k=1 k=1 
or 
Fy iy = — Md? y/dt?, (23) 


for the values of FR and & chosen. 

For vector spatial ballistocardiography we use the position vec- 
tor r and consider the spatial forces to be denoted by F. The sub- 
scripts are given the same meaning as in the one dimensional 
case. Instead of equations (1) and (2), we have 


d(m, dt, /dt)/dt = F,, + F,,, (24) 
and 
d(m,dt,/dt)/dt=F,,+F,,. (25) 
Instead of equation (28), we have 
Fe Md?r,/dt?. (26) 


Discussion. The suspended table type of low frequency ballis- 
tocardiograph (Burger et al., loc. cit.) or the aperiodic bed mercury 
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ballistocardiograph (Talbot et al., 1954; Deuchar et al., 1955) 
would appear to furnish desirable characteristics. In the low fre- 
quency ballistocardiograph & is relatively small and & can be made 
as large or almost as small as desired. In the aperiodic bed & is 
zero, but R has a definite magnitude. Both types of machine have 
been shown to give similar records. 

It is of interest to examine for a given set of conditions the 
variation of 1/\/R2/@2 M2 + (k/w2M ~1)2 and of arc tan whe/(k - Mw”) 
for different values of w. We do not know the value of F for a 
mercury table. However, & being equal to zero would tend to make 
are tan wR/(k - Mw”) smaller, but would tend to make 
1/VR?2/a2M2 + (k/w? M — 1)? closer to unity. 

For the low frequency, ballistocardiograph representative values 
of the constants are 


M =8 x 10* gm, 

k = 3x 104 gm/sec?, 
and 

R = 1.8 x 10* gm/sec. 


The length of the supporting cables is assumed to be about 250 
cm in calculating & and there is assumed to be about 85% over- 
shoot in calculating R. The lowest frequency component of any 
force produced by the heart is about 1 cycle/sec (heart rate of 60 
beats per minute). Hence for some purposes we need not concern 
ourselves with frequency components lower than 1 cycle/sec, par- 
ticularly if respiration is suspended. Under these conditions the 
maximum phase shift is represented by 0.006 seconds of time. The 
value of 1/\/ R2/w2 M? + (k/w? M — 1)? fora frequency of 1cycle/sec 
ARES For frequencies greater than 1 cycle/sec, 
1/VR2/@2 M2 + (k/w? M — 1)? is even closer to unity. 


Summary. 1. Ballistocardiography may be considered as an ex- 
ample of the two particle problem in physics where the particles 
are of variable mass. 

9. The internal forces are those produced by cardiac activity and 
exerted by the blood and cardiovascular structures on the remainder 
of the body. 

3. Physical analysis shows that the internal forces not only 
arise from moving blood but are a resultant of the motions of blood, 
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cardiac muscle, vascular walls, pulmonary tissue, and all other 
structures moving relative to the rest of the body or a reference 
system based on the skeleton’s position. 

4, The external forces are those exerted by the table supporting 
the patient’s body. 

5. The basic problem in ballistocardiography is the determina- 
tion of the magnitude of the internal forces as a function of time. 

6. The suspended low frequency table and the mercury bed are 
promising instruments. 

7. The acceleration curve of the low frequency table is linearly 
related to the various components of the internal force irrespective 
of the frequency components of the latter. 


The authors wish to thank Professor N. Rashevsky and Dr. 
George Karreman for critically reading the manuscript and offering 
valuable suggestions. 
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The probabilistic theory of random and biased nets is further developed 
by the ‘‘tracing’’? method treated previously. A number of biases expected 
to be operating in nets, particularly in sociograms, is described. Dis- 
tribution of closed chain lengths is derived for random nets and for nets 
with a simple ‘“‘reflexive’’ bias. The ‘tisland model’’ bias is treated for 
the case of two islands and a single axon tracing, resulting in a pair of 
linear difference equations with two indices. The reflexive bias is ex- 
tended to multiple-axon tracing by an approximate method resulting in a 
modification of the random net recursion formula. Results previously ob- 
tained are compared with empirical findings and attempts are made to ac- 
count for observed discrepancies. 


There is increasing interest in the theory of linear graphs which 
is felt to be a promising mathematical framework for certain types 
of biological and sociological theories, in particular those dealing 
with structure. Indeed a graph is primarily a representation of 
structure in its most essential sense: the points or nodes are the 
elements devoid of all content; the lines connecting them represent 
relations which in the simplest case can be characterized by either 
presence or absence between each pair of elements. On the simplest 
level there is an isomorphism between a linear graph and a matrix. 
Also the operation of matrix multiplication is applicable to the de- 
duction of certain theorems about graphs. This mathematical ap- 
paratus has been utilized to some degree by several authors (Landahl, 
1947; Landahl and Runge, 1946; Luce, 1950; Luce and Perry, 
1949; Shimbel, 1948). 

Topological properties of linear graphs have found application 
to problems of general biology dealing with the ‘‘flow chart’’ or- 
ganization of physiological function in the work of N. Rashevsky 
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(1954, 1955a,b, 1956) and E. Trucco (1956). Applications to social 
psychology have also been proposed (Harary and Norman, 1953; 
Cartwright and Harary, 1956). 

In all of this work, the theory of linear graphs appears as an 
‘‘exact’’ theory in the sense that each graph is supposed to de- 
scribe the complete structure of some entity at least in its relevant 
aspects. In dealing with structures of aggregates involving very 
large numbers of elements, it seems advisable to develop statisti- 
cal, probabilistic, or stochastic theories. Some beginnings in this 
direction have been made (Solomonoff and Rapoport, 1951; Rapo- 
port, 1951a,b, 1953a,b; Katz and Powell, 1954; Allanson, 1956). 

However, even a probabilistic theory can be either ‘‘exact’’ or 
approximate. In an ‘‘exact’’ probabilistic theory it is necessary 
to define rigorously a ‘‘sample space,”’ i.e., the set of all elemen- 
tary events with their a priori assigned probabilities, of which the 
events of interest (where probabilities are to be calculated) are 
compounded. Such an approach is taken, for example, by L. Katz 
and J. H. Powell (loc. cit.). The fundamental event is the exist- 
ence or non-existence of a relation between each ordered pair of 
elements in a set of N elements. Since there are N(N —1) such 
pairs from which set 2%°¥—» subsets can be selected, it follows 
that Katz’s sample space is the set of 2’“%- linear graphs. If 
the probability for the existence of the relation between each pair 
of the set of N elements is fixed, a set of probabilities is induced 
on each set of points of the sample space, from which the statisti- 
cal parameters of such linear graphs can be calculated. Restric- 
tions, such as fixing the total number of existing relations, lead to 
the consideration of a subspace of the total sample space, which 
is in fact considered by Katz and Powell. 

We could also begin with another assumption, for example, that 
exactly a directed line segments (‘‘axons’’ in our terminology) 


issue from each of N elements. Since there are (M4) ways of select- 


ing the targets of each set of @ axons, it follows that ban graphs 


of this sort are constructible, indeed all equiprobable if every ele- 
ment is an equiprobable target of every axon. In both of these ex- 
amples, it is tacitly assumed that all elements are distinguishable. 
If they are not, the probability distributions on the sample spaces 
must be modified accordingly. This straightforward ‘‘exact’’ ap- 
proach depends on combinatorial analysis for its mathematical 
machinery and soon leads to formidable difficulties. 
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Another approach is to start not with the sample space of all 
possible graphs but with the probabilities of the events in the con- 
struction of a graph. Starting with an arbitrary element, there will 
be a certain probability distribution for all possible targets of its 
axons and so on for each target thus reached. If we consider this 
‘“‘tracing’’ of the graph as a stochastic branching process, we are 
still within the realm of an ‘‘exact’’ probabilistic theory. It is 
possible, however, to abandon rigor in the interest of a more man- 
ageable model by dealing only with the expected events at each 
step of the tracing instead of with the entire distribution of all 
possible events. This was the method adopted by the author in 
previous papers on random and biased nets (Rapoport, 1951, 1953a,b). 
The relation of this method to the ‘texact’’ stochastic process was 
discussed by H. G. Landau (1952). In the special case where only 
a single axon issues from each element, the ‘‘exact’’ theory can be 
easily handled. In previous papers (e.g., Rapoport, 1948; Shimbel, 
loc. cit.) the distribution of cycles of all lengths in a random net 
have been calculated. 

We begin our present investigation by calculating the distribu- 
tion of the lengths of ‘‘closed chains’’ in single axon random nets. 


Closed chains in single axon random nets. The construction of 
such a net begins with the arbitrary selection of an element and 
the tracing of its axon to its target, say, by a random device in 
which all elements are equiprobably represented. The tracing then 
continues to the next target and so on until a previously contacted 
element is contacted for the second time. When this happens the 
tracing ends and the chain is closed. Required is to find the 
probability distribution of closed chain lengths. 

Let p(n) = the probability that the tracing is not yet ended on 
the nth step, i.e., that n+ 1 distinct sites have been occupied. 
Then if the total number of sites is N, we have 


pin) = 28, (1) 


where N, = N(N- 1)...(N-n+1). Then the probability of a 
chain of length exactly n must be 


N N, N N-n nN, 
pelea eres PE Tela N )< yar ® 
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Remark. It should be noted that while this result can be ex- 
pected to reflect the distribution of chain lengths in a large popu- 
lation of nets constructed in this manner, it cannot be expected to 
reflect the distribution of chain lengths in a single net, because of 
the overlap of chains. Thus the existence of a closed chain of n 
steps implies the existence of a closed chain of n — 1 steps (namely, 
the one which starts with the second step of the first chain), etc. 
In particular, if the entire population happened to fall apart into a 
number of such separate chains of equal length, the distribution of 
chain lengths in that population would be an equiprobable one 
(‘‘square’’) and not given by (2). This inapplicability of the dis- 
tribution of chain lengths to a particular population in which an 
axon issues from each site is important to keep in mind when we 
try to apply such theories to actual biological or sociological nets. 

Consider now a sociogram in which every member of a population 
names one ‘‘best friend.’? The result is a single axon graph of the 
sort described. However, even if the ‘‘best friends’’ are named 
entirely at random, the chain lengths in any particular population 
cannot be expected to be given by (2), because of the dependencies 
mentioned. Equation (2) describes the distribution of chain lengths 
where one chain is arbitrarily selected from one population in a 
large set of populations. 

The independent structures in a Single random net of this kind 
are the ‘‘cliques.’’ Consider a net in which exactly one axon is- 
sues from each site. Therefore each axon has a unique target. 
However, more than one axon may reach the same target. Some 
sites may not be targets of any axon. Thus if we trace the chain 
“‘forward,’’ that is, to successive targets, we get no branching and 
the chain ends when a site is repeated. If we trace ‘‘backwards,”’ 
that is, to the origins of the axons, we may get branchings but no 
repeated sites. Each branch will end when an ‘‘isolate’’ is reached, 
i.e., a site which is not a target. Thus if we perform both a for- 
ward and a backward tracing, we reach a certain subset of the 
population with which no connections exist from other subsets 
either via forward or via backward tracings. We call such subsets 
“‘cliques.’’ Figure 1 shows such a clique for a single axon net. 

Thus cycles, closed chains, and cliques are structures charac- 
teristic of single axon nets. A cycle is completely described by 
its length (except for the identification of the sites involved). A 
chain consists of a cycle and a tail and is completely described 
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by its total length and the length of its cycle. The problem of 
describing a clique is more complicated and bears an obvious re- 
lation to similar structures in organic chemistry (cf. Pélya, 1936). 

Once cliques are classified in some way, the problem arises 
concerning the distribution of the various types of cliques in a 
single axon random net. These distributions should reflect the 
actual frequencies of occurrence (provided the number of cliques 
is not too small) since the cliques are independent of each other. 

In multiple axon nets the substructures become increasingly 
complex, and we are forced to concentrate on simpler character- 
istics of such nets. One approach taken by the author (Rapoport, 
1951a) was to calculate the recursion formula for the expected 
number of new sites contacted on each step of a tracing, that is, 
the expected number for each step is calculated on the assumption 
that the expected number was actually contacted on the previous 
step. In another paper the actual probability distribution for the 
number of sites was calculated for a given step, given that the 
total number of sites previously contacted was known (Rapoport, 
1951b). Elsewhere the expected total number of sites contacted 
during the whole tracing, expressed as a fraction of an infinite 
population, was derived as a function of the ‘‘axon density’’ (num- 
ber of axons issuing from each site) (Solomonoff and Rapoport, 
loc. cit.). 

All these results were obtained for random nets where the under- 
lying assumption is the equiprobability of targets for each axon. 

In a situation where we may hope to apply the theory (e.g. struc- 
ture of neural nets, sociometry, epidemiology, chain reactions) we 
should expect the operation of substantial biases: the targets of 
axons are not equiprobable. The introduction of bias leads to 
formidable difficulties because the abandonment of the assumption 
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of equiprobability destroys the independence of certain fundamental 
events. In particular the steps of the tracing cease to be inde- 
pendent. Nor can they in general be reduced to a Markoff process 
except under special assumptions, for now the whole history of a 
tracing has a bearing on the probability distribution of the events 
associated with each step. 

However, besides the difficulties of calculation, there are other 
difficulties involved of a more fundamental nature, namely, of 
stating precisely in mathematical terms what bias is operating. 
We wish naturally to investigate biases which can be reasonably 
expected to reflect some situation of interest, but then we find 
that it is difficult to translate intuitively conceived biases into 
formal postulates. 

To take an example, it is a commonplace observation. that the 
relation of acquaintance is to a certain degree transitive in the 
sense that acquaintances of acquaintances have a greater proba- 
bility of being acquaintances than arbitrarily selected individuals. 
However, this ‘‘bias,’’ although seemingly clear, is not at all easy 
to express in terms of biases imposed on the events associated 
with the tracing procedure. Therefore, it is advisable to introduce 
biases formally as affecting the events associated with the tracing 
procedure, so that they can enter the calculations, instead of try- 
ing to derive them from the actual situation envisaged. We hope 
that some of such formally introduced biases may be approximately 
equivalent in their effects upon the derived statistical parameters 
to the biases actually operating. 

We submit a tentative list of such formal biases. 

1. Geographical distance bias. The elements are supposed to be 
immersed in a metric space, and the probability that an axon is- 
suing from one element ‘‘synapses’’ on another is some function of 
the distance between the two elements. Aside from the difficulties 
concomitant to this bias in the calculation of conditional proba- 
bilities, this model poses the problem of choosing a suitable ‘‘space”’ 
for a given situation. Clearly, if the net is a sociogram, we are 
dealing with ‘‘social space.’’ Distance in such spaces could be 
appropriately defined, but nothing whatsoever is known about the 
dimensionality or metrics of spaces which are appropriate models 
of social space. In principle we could pose the problem of in- 
ferring the character of such a space from the observed statistical 
parameters, but a theoretical framework for this approach is still 
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lacking. Where ordinary geometric space can be taken, the latter 
difficulty is obviated, but the former one remains, as is seen in 
some formulations of mathematical genetics of populations. 

2. The ‘island model.’? We may suppose that the elements are 
divided into a number of mutually exclusive subsets within which 
the connection probabilities are random but where a finite proba- 
bility exists of crossing from one such subset to another. This is 
the ‘island’? or ‘‘cellulated’? model which has been treated in 
some detail in mathematical genetics (cf., e.g., Wright, 1932). It 
seems a reasonable model for some types of sociograms. 

3. Overlapping acquaintance circles. Here an element is generally 
a member of more than one set, so that even if no crossing of set 
boundaries occurs, that is, no two elements are connected unless 
they belong to the same set, nevertheless the entire population can 
be connected. This model underlies some of the author’s previous 
results on biased nets (Rapoport 1953a,b) in which very crude ap- 
proximations were used. The following problem illustrates this 
type of model in a simple but interesting case. 

Let the elements be arranged in a linear array and each axon 
have equal probabilities of connection to each of its 24 neighbors 
(k on each side). It is required to compute the distribution of the 
lengths of closed chains and cycles. Put in another way, consider 
a random walk in which 1 to & steps may be taken in either direc- 
tion and which terminates once a site previously occupied is oc- 
cupied again. Compute the distribution of the lengths of walks. 
For k = 1, the problem is trivial, because the walk ends if and only 
if the direction is reversed. For & = 2, the problem is considerably 
more involved. The solution for this case was obtained by C. 
Foster and the author (in press). 

The solution is obtained in iterated form as a sequence of ‘‘con- 
volutions.’’ Define 


i=[x/3] ; (x - 2%)! 
F, (@) = me (i/2)"t* Toe (3) 
Then if 
AAD = », F, (x -#) (1/4)! + F,(2), (4) 


i=2 
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then the probability of a path of at least n steps is given by 


p(n) =2F,(n), (5) 
and the probability of a path of exactly n steps by 
P(n) = 2[F ,(n - 1) - F,(n)I. (6) 


Generalization is indicated to & > 2 and to grids of several di- 
mensions. If solutions for any & and any number of dimensions 
were available, we could compare such distributions with actual 
ones (say, in sociograms where the ‘‘step’’ is the naming of the 
best friend). The best fitting model, where & and the number of 
dimensions would be the parameters, could then be taken as an 
‘‘overlapping acquaintance circle’? model of social space, an ac- 
quaintance circle being a neighborhood of a certain size in a space 
of a certain number of dimensions. 


4. Reflewive bias. This is a bias of connection from a target 
back to the origin of an axon. This model also seems appropriate 
in the tracing of a sociogram. More generally, the reflexive bias 
operates with respect to all elements already reached in a tracing, 
the magnitude of the bias being a function of the number of steps 
to the element in question. 


5. Force field bias. Here certain elements are endowed with a 
‘‘drawing power,’’ so that connection to those elements are more 
probable. In a special case we may suppose that the population 
consists of two classes, the ‘‘popular’’ and the ‘‘unpopular,’’ the 
popular elements having a greater probability of being targets. In 
a neural net these elements would constitute some ‘‘center’’ to 
which most connections go. 

In this paper, we consider the following models related to the 
probabilistic characteristics of biased nets. 

a) Distribution of chain lengths in single axon tracings with 
mutually exclusive sets with a finite probability of crossing over 
for the case of two sets. 

b) Distribution of chain lengths in a single axon tracing with 
reflexive bias. 

c) An approximate recursion formula for the number of elements 


reached in successive steps of a multiple axon tracing with re- 
flexive bias. 


=e 
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d) A recursion formula for the case of a force field bias for the 
special case of two classes. 

Data from a sociogram will be compared with the theoretical re- 
sults previously obtained and discussed in the light of the theoretical 
results obtained in c and d. 


Single axon tracing in a two-set population. Let the two equal 
sub-populations be called A and B. Let the population of each be 
M. Let p(n, &) be tke probability of tracing a chain of at least 
n — 1 steps, i.e., one in which n sites are occupied, of which k Sn 
are in A, hence n— are in B. Let the probability of crossing from 
A to B or vice versa be xz. At the completion of the (n — 1)st step, 
the last site occupied is either in A or in B. Let the correspond- 
ing probabilities be p, (n, &) and pz (n, k) so that 


p, (nk) + Pp(n, k) = p(n’). (7) 


Clearly, then, the probability of getting a chain of at least n- 1 
steps will be given by 


p(n) = )° p(n), (8) 


k=0 


and it remains to calculate the p(n,k) which in turn will be ob- 
tained if p, (n,*) and p,(n,*) can be determined. These latter pair 
of functions of two discrete variables satisfy the following pair of 
difference equations 


M-(k-1 M-(k-1 
Bee eie t yiha-Do 
M -(n —k) M-(n-k) (2) 

Da(n+1,k)=py(n,) & Sl ypsto + pg(n,k) (1-2) ——— 


We will justify the first of these. The left side is the probability 
that (n + 1)st site to be occupied is in A, there being & occupied 
sites in A. This can happen if the nth site was also in A, in 
which case there were & — 1 occupied sites in A on the preceding 
step. The probability of this happening is compounded of 
p,(n,k — 1), A - z) which is the probability of non-crossing, and 
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[M —(k —1)l/M which is the probability of not terminating the 
chain. Hence the first term on the right. Otherwise, the nth site 
was in B, in which case the path crossed from B to A; hence the 
second term on the right. The second equation is justified by a 
similar argument. Adjoining to the system (9), we have the initial 
conditions 


Pa (1, 1) = 1/2; Pa (1, 0) = 0; 
pp (1,1) = 07 p, (1,0) ="172. 


(10) 


Equations (9) and (10) determine the functions p, (n,%) and pp(n,k), 
hence p(n, k), hence p(n), which for the special case = 1-2 =1/2 
reduces to the case of the random single-axon net treated previ- 
ously (Rapoport, 1948). The case a = 1 is equivalent to that treated 
by Shimbel (loc. cit.). The latter is a model for the case where 
choices are arbitrary but always, say, of the opposite sex. Gen- 
eralization to an arbitrary number of sets leads to a much more 
complicated system in which the single parameter & must be re- 
placed by a set of integers representing a decomposition of n. 
The solutions of the system are then functions of n and m-1 
additional parameters where m is the number of sets. To obtain 
p(n), these functions must be summed over all the decompositions 
of n. 


Single axon tracing with reflexive bias. In a random net the 
probability that one of the n sites already occupied is a target of 
the axon issuing from the last site occupied is 


n 1 1 1 

VON ee oe (n terms). (11) 
The equality of the terms on the right constitutes the randomness 
condition. If we weight these probabilities with factors sere 


(¢=1, 2,...,n), we have reflective bias. We must, of course, 
have 


N 


Dd, By-i= (12) 


ial 


and it is most natural to suppose that 8 ,,_; 18 independent of n and 
is a monotone decreasing function of (n-2). Call it B(i). A 


~~ one 
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straightforward generalization of (1) then gives 


WV BEEN =BC0)= AOI. ye T 
t=1 


p(n) = - 


N"™ 
k 
De B (i) 


n i=1 


J, N 


(13) 


For large N, we may suppose that N is infinite but that B(i)/N = 
5(z) is finite so that 


) ‘gale (14) 


Then 


n k 4 
log p(n) = es log . oe val : (15) 


k=1 t=1 


or, passing to approximations by integrals, 


lop p(n) -{- lap f -[ ba) de es (16) 
i?) 0 


Suppose now 4() is a decaying exponential, so that by our normal- 
ization condition (14), we must have 


h(x) =ae **. (17) 
Equation (16) becomes 


n k n 
log p(n) -{ log A = | ates | dk ef log e~°* dk 
0 0 0 


Eid Sepa es 20718 
0 


(18) 


so that 
Ray eee (19) 
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In other words, in the special case where the reflexive bias is an 
exponential function of the ‘‘distance traversed,’’ the probability 
of having a chain of length at least n falls off with n as the posi- 
tive half of the gaussian curve, whose variance is the reciprocal 
of the decay constant of the distance function. 


Reflexive bias in multiple axon tracings. Formally, the result 
just obtained can be applied to a multiple axon tracing, but its 
interpretation in that case is far from clear. In the single axon 
tracing each axon has a unique origin; thus in a chain traced for- 
ward each site reached after the first has a unique antecedent. In 
a multiple axon tracing the sites reached on the nth step are clearly 
antecedents to those reached on the (n + 1)st. But it is not at all 
clear which are antecedent to which. This information is lost 
where the tracing is described only in terms of the numbers of 
sites (or fractions in an infinite population) reached at each step. 
Nevertheless we will extend the results of the previous section to 
multiple tracings. We recall the recursion formula derived for a 
random net (Rapoport, 1951a), 


t 
P(¢+1)= |1- rh P(j)| (1 = e~ 2? OF, (20) 
j=0 
where P(t) is the fraction of sites newly contacted on the ¢th step, 
and a is the number of axons issuing from each site. In applying 
this formula, we can solve for a, which appears formally as a func- 
tion of ¢, thus 


t 


i ead 


1 j=0 
a(t) = P® log esa (21) 
1-1 PG) 
j=0 


To the extent that this function of ¢ appears to be a constant, the 
randomness hypothesis for a given net which yields the values 
P(j) in a tracing is corroborated. It was further shown (Rapoport 
1952a) that the bias of overlapping acquaintance circles is ex- 
pected to yield 


a(0) =a, 
a(t) & a’ < a(t >.0). (22) 


So es 
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That is, «(¢) drops from a to a lower value for ¢ = 1 as calculated 
from P(1) and P(2) and remains approximately at that level for the 
successive steps. Certain data from an information spread ex- 
periment (Rapoport, 1953a,b) indicated that a«(1) is indeed drasti- 
cally less than (0) in all cases but that subseauently a(¢) tends 
to increase with ¢. This increase was attributed to the conditions 
of the experiment and was expected to disappear in a sociogram 
tracing if the overlapping set model is valid. We now examine the 
behavior of a(t) in the reflexive bias model. 

In a random net the probability that a given axon finds a target 
in any subset is proportional to the size of the subset. In particu- 
lar the probability that this target is in P,, the fraction newly con- 
tacted on the ith step, is P;. Let now a bias be imposed so that 
the probability of finding a contact ingests Pak let B,_,)» Where 
Boe. depends on the number of steps between the ¢th step (the step 
in question) and the zth step. Then a bias y is imposed on the re- 
maining targets, i.e., those sites not yet contacted. The normali- 
zation condition demands that 


Fare. 0-2 Pe (23) 


i=0 i=0 


Solving for y, we have 


iL y Pies) 


oo 


i=0 


y= (24) 


We note that for B,_,=9, y reduces to unity, as should be the 
case. Recall now that equation (20) is an exponential approxima- 
tion to 


t t fe NaP,(1 = 2.)] 
resn-(-2F, 1-|1-N be) fe (25) 
. in ' 


5-2 


Here the second term in the brace is the probability of noz being a 
target of any of the axons going to unoccupied sites. When the 
bias y is imposed on these axons, the exponent of the bracket ac- 
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quires y as an additional factor. Passing again to the exponential 
approximation, we have for this case 


ACs UNE MaeDPIA NERA FO) (26) 


j=0 


Remembering that o(¢) is the coefficient of P, in the exponent, and 
substituting (24) into (26), we obtain 


YB; 
t=O 


ath of Is ae ee (27) 


t 


Eee: 


i=0 


The denominator in parentheses always decreases. The numera- 
tor, on the other hand, may increase or decrease depending on the 
behavior of P, and on the bias parameters 8,_, chosen. It is clear 
however that, for B,_; > 0, ¢ > 0, a(t) < a, as is corroborated by all 
data examined so far. We may ask under what conditions will «(¢) 
be constant for ¢ > 0, so that both the overlapping set model de- 
veloped previously (Rapoport 1953b) and the present reflexive bias 
model are indistinguishable. This condition is obviously 


nA FB Al Y ") (28) 


i=0 i=0 


where A is a constant, a system of linear equations which may be 
solved for B,_;. We would expect B,.; to be a monotone decreas- 
ing function of 7 and it would be interesting to verify this result in 
these cases where &(¢) is approximately constant for ¢ > 0. 


The ‘‘force field bias.’’ This case is actually a generalization 
of the ‘‘mutually exclusive set’? case where now each set or ‘‘is- 
land’’ is characterized by a different probability of being crossed 
into. In a still more general case, we have a set of directed cross- 
ing probabilities between each ordered pair of sets. The simplest 
case is where a ‘‘center’’ exists, such that all axons find targets 
in the center. That is to say, the probability of crossing into the 
center is unity; the probability of crossing out of it is zero; and 
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all targets are equiprobable within the center. Except for the ar- 
bitrarily selected set Py, all the P, will be in the center; and from 
the first step of the tracing on, the model is equivalent: to a ran- 
dom net. However when dealing with actual data the existence of 
the center is not known. It may be detected by observing the be- 
havior of «(¢) which is deducible directly from the P,. This ‘‘ap- 
parent’’ «(¢) will no longer be constant. Another function will be 
constant, which we now compute. 

Let f be the fraction of elements in the ‘‘center.’’ Then, since 
P, is chosen arbitrarily, fP, of these are expected to be in the 
center. All the subsequent P, are also in the center. Equation 
(20) now has the form 


P(t +1)=\|f-fP, - ae Py ties od: (29) 


bee 


Solving for the coefficient of P, in the exponent as before, we have 


us P; 
pee te SS (30) 
a = — log —————_ 
: i 2 ttl 3 
f, os a P, 
i=0 
where now 
freafie tt =). (31) 


We see that f, < 1 so long as f<1. To compute theoretically the 
behavior of the ‘‘apparent’’ o(¢) given by (21), the P, must be de- 
termined from the recursion formula (29), a laborious task since no 
closed form is known for P(j). However, the effect of the ‘‘cen- 
ter’? can be estimated in any given set of data in which the P(j) 


are known. 


Applications of the foregoing to some experimental data. A so- 
ciogram of an elementary school of some 250 children was obtained 
in the following way. Each pupil was asked to name ten friends, 
who were also pupils in the school, in the order of intimacy. In 
this way there was obtained a card for each child on which ten 
names were listed, the first of which was supposedly the “‘*best 
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friend” and so on. Single axon and multiple axon tracings could 
then be made from these cards. We could, for example, make sev- 
eral single axon tracings through the ‘‘best friend’’ in which the 
target of each axon was always the first name on the card, or 
through the Ath friend, where the target was always the Ath name. 

So far the only multiple tracings were made with a =2. In par- 
ticular it had been conjectured on the basis of the overlapping set 
model that the behavior of &(¢) should be given by equation (22). 
Moreover, it was argued that tracings through the names high on 
the list should result in a greater departure of a(t) from a for ¢ > 0 
than tracings through names low on the list. Accordingly 10 trac- 
ings each were made through the first two names, the middle two 
names, and the last two names on the list. The P(t) averaged over 
the ten tracings and the computed a(¢) are given in Tables I, II, 
and III. 


TABLE I 


Tracing through the first two names on each card. 


TABLE I 


Tracing through the middle two names on each card. 


TABLE OI 


Tracing through the last two names on each card. 


We note the following characteristics of the data: 


1. «(0) = 2 in all cases, agreeing exactly with the fixed value 
of the parameter a, which should be the case since P(0) are se- 
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lected at random. This is reflected in the near-equality of P(1) in 
all three cases. 

2. «(1) <2 in all cases; moreover the drop is most pronounced 
in the tracing through the first two names, less pronounced in the 
tracing through the middle two names, and least pronounced in the 
tracing through the last two names on the list in qualitative agree- 
ment with the theory. 

3. Instead of the steady rise in a(t) for ¢ > 1 which was observed 
in the information spread experiment (cf. Rapoport 1953b), «(¢) de- 
clines in the first half of the process and rises in the last half. 
(Except for the last value in Table I.) 

4. We note that P(t) increases in the first part of the process 


t 
and decreases in the second. The quantity log ( =i PU 


tt1 
( - a Pea| , whose calculation is not shown in the tables, be- 


haves similarly. By our theory, the ratio of these two quantities, 
that is, a(t), should remain approximately constant fort > 1. Their 
similar behavior facilitates this result, but apparently the varia- 
tion of P(¢) more than compensates for the variation of the other. 

If the concave upward curve for a(t) is real and not an artifact 
of random fluctuation, then the overlapping set model is not cor- 
roborated by these data, and we have the choice of either seeking 
another model or of trying to modify the existing one by some addi- 
tional considerations. Two attempts were made in this direction, 
one involving a reflexive bias to the immediately preceeding step 
only and one involving the postulate of a ‘“center.’’ Both modifi- 
cations accounted for the successive drops in a(t) for ¢ = Core: 
but resulted in further drops for the remainder of the process thus 
failing to explain the subsequent rise. 

Needless to say additional modifications of the model become 
increasingly ad hoc and eventually lose theoretical significance. 
Modifications are worth while only if on their basis additional 
effects can be predicted which are then verified. 

It is possible to account qualitatively for the initial fall and 
subsequent rise of a(t) on the basis of an ‘tisland’’? model and to 
derive an additional prediction to be tested. Suppose the popula- 
tion is divided into several mutually exclusive sets, and suppose 
P(0) is sufficiently small so that not all sets are represented in 
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that sample. While no new sets are represented in the P(t), the 
tracing is equivalent to one in a random net but with a smaller 
population. As a result, the apparent o(¢) will be declining. This 
has not been shown theoretically although it is intuitively evident 
and can be demonstrated on our set of data. When, however, at 
some stage of the process a new set is ‘‘seeded,”’ the initial rate 
of increase of P(t) in that set is great, as it always is in the be- 
ginning of the tracing and so the over-all rate of increase of P(¢) 
gets a ‘‘boost’’ (or its rate of fall declines) with the result that 
the apparent o(¢) again increases. This should occur every time a 
new set is crossed into while old ones become saturated. If this 
is so, then the fluctuations in «(¢) should be the more pronounced 
the smaller the initial fraction P (0). 

To test this hypothesis, tracings were done through the first two 
names, with P(0) =.04 instead of .02 as in the other tracings. 
The results averaged over five tracings as shown in Table IV. 


TABLE IV 


Tracing through the first two names with a larger initial sample. 


The ‘‘dip’’ in the middle is still present though seemingly less 
pronounced. The last value is again lower than the preceding one 
as in Table I, with which Table IV is in fairly good agreement. 

The question of the ‘“‘dip,’’ then, remains open. Without a 
thorough statistical analysis it is impossible to say whether the 
effect is real or a fortuitous result of fluctuations. Thorough sta- 
tistical analysis cannot be undertaken without a fairly rigorous 
mathematical model of the process, which, in turn, has to be a 
genuine stochastic model rather than a recursion eauation for the 
most likely values of P(t). ; 


Discussion. It is evident that questions relating to the sta- 
tistical structure of nets are exceedingly involved both in their 
formulation and in the computations they entail. In particular all 
the multiple axon tracing models so far considered are devoid of 
mathematical rigor even in their simplest formulation and become 
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increasingly more crude as they are modified by imposed biases, 
so that it is impossible to tell at which point a mathematical theory 
gives way to a weakly rationalized, ‘‘guessed’’ equation which 
might as well have been obtained empirically. 

On the positive side, a large number of statistical parameters 
can be observed in any real net. Given a certain stability of these 
parameters with respect to several nets of the same type, we could 
hope to pass from an empirical observation of their behavior to the 
building of mathematical models from which the statistical param- 
eters could be derived. 

As an example, consider once more a single axon sociogram. If 
it is a random net, then the distribution of sites having respec- 
tively none, one, two, etc., axons terminating on them ought to be 
a Poisson distribution with unity as mean, i1.e., 


p(n) =e “/nt, (32) 


where p,(n) is the probability of a site receiving n axons. Sim- 
ilarly in a tracing with a axons per site, the distribution ought to 
be 


p(n) =a" eo SM (33) 


We have examined several single axon tracings in our school 
sociogram. The single axon tracing through the ‘‘best friend’’ 
gave a very good fit to the Poisson ‘distribution. However all 
other tracings (through ‘‘second best friend,’’ etc.) gave depar- 
tures from the Poisson distribution, all in the same direction—the 
number of ‘‘normals’’ (those receiving a single choice in a single 
axon tracing) was reduced in each case. The ten-axon tracing in- 
volving all the friends named gave a definitely bimodal distribu- 
tion with a depression at the expected peak. Here may be evi- 
dence of a bias roughly described by the force-field model: the 
population can be divided into the popular and the unpopular. It 
should be noted incidentally that the Poisson distribution charac- 
teristic of a random net is sensitive to some biases but not to 
others. It would not, for example, be sensitive to the bias of the 
‘island model’’ provided the populations of the islands were suf- 
ficiently large. On the other hand a departure from the Poisson 
distribution should become apparent in a “force field’? bias model. 

Another striking result of the single axon tracings has to do 
with the number of resulting cliques and the sizes of their cycles. 
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The ‘“‘best friend’’ tracing resulted in our population of 250 falling 
apart into 43 cliques, 12 of them couples, all cycles but one being 
2-cycles. Tracing through the ‘‘second friend’’ gave only 27 cliques; 
tracing through the last friend named only 12 cliques, cycles of 
various sizes being well represented. Theoretical explanation of 
this decrease of the number of cliques with the ‘‘social distance”’ 
associated with the axon which is traced awaits the creation of an 
appropriate mathematical model. 

A method is thus indicated for describing nets of relations in 
terms of their departures from random nets. Different statistical 
parameters show different ‘‘sensitivities’’ to departures resulting 
from different kinds of bias. It seems desirable to have an ex- 
tensive theory of these dependencies so that a whole battery of 
measurements can be applied to a given net to determine what sort 
of bias predicts the greatest number of relations among its sta- 
tistical parameters, expecially since a great number of such sta- 
tistical parameters can now be empirically obtained by proper use 
of computing machinery. 
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THE USE OF ORGANIC DYES IN ELECTRON MICROSCOPY 
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The extinctions produced by several organic dyes in the light micro- 
scope and the electron microscope are compared. The two have the same 
order of magnitude. The influence of the methacrylate embedding medium 
is not great enough to hamper visibility in the electron microscope. It is 
pointed out that a wide variety of organic stains should be useful for 
electron microscopic studies. 


Osmium tetroxide has been the fixative and stain of choice for 
biological specimens for electron microscopy since 1953. How- 
ever, it has long been known that the use of osmium limits the sub- 
sequent application of dyes and enzymes in light microscopy both 
for strictly morphological and for cytochemical purposes (Romeis, 
1948, p. 66). Consequently, for the most part, only indirect means 
can be used for the localization at submicroscopic levels of baso- 
philia, DNA, enzymes, glycogen, etc., when osmium fixation has 
been used. 

When tissues are prepared for electron microscopy by freezing 
and drying (Gersh et al., 1957), they are stainable by a variety of 
procedures. The results of staining for glycogen (Bondareff, 1957) 
and for RNA and DNA by gallocyanin chromalum or by the Feulgen 
stain (Finck, 1957), which give characteristic appearances in the 
light and electron microscope, make it clear that descriptions of 
localization may be made at the submicroscopic level. The work 
of the foregoing authors shows that the stains employed had suffi- 
cient contrast for electron microscopic observation. It therefore 
becomes important to examine from a theoretical viewpoint the 
problem of contrast changes when organic substances are used as 
electron microscopic stains. 


*Present address: Institute for Muscle Research, Marine Biological 
Laboratory, Woods Hole, Massachusetts. 
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This paper will not deal with general problems of cytochemistry 
in electron microscopy. Nor will this paper attempt to review or 
extend the principles involved in the production of images in an 
electron microscope. Although peripheral questions will be dis- 
cussed, this paper will attempt only one major theoretical thesis, 
namely, that common histological organic dyes may profitably be 
used as electron microscope stains. When so used sections of 
stained tissue may be compared with unstained controls and differ- 
ences in density will be apparent either on the electron microscope 
screen or in electron micrographs. By such comparisons we may 
attempt to draw conclusions, as did Bondareff (loc. czt.) and Finck 
(Loc. cit.),as to the chemical composition of those tissue elements 
which the dye has stained. 

For the development of such a thesis the simplest possible case 
will be considered. It will be assumed that a sufficiently thin sec- 
tion is being examined in an electron microscope with a suffi- 
ciently small objective aperture so that the intensity dimunition of 
the electron beam follows a simple exponential law. Phase con- 
trast effects will be ignored, as will contrast changes due to the 
use of specific high or low contrast film in taking electron micro- 
graphs. 

In considering contrast in the light microscope, only light ab- 
sorption will be considered to be significant. 

Since nothing is known concerning the quantitative submicro- 
scopic distribution of protein, nucleoprotein, or dyes, compara- 
tively little can be achieved by a direct calculation of contrast in 
the electron microscope. The problem is one of not knowing the 
concentration or density of the tissue elements or stained tissue 
elements under consideration. This problem will be circumvented 
by comparing the extinction of tissue similarly stained for optical 
(light) viewing and for electron optical viewing. Such a compari- 
son will also permit us later to refer the criterion for visibility in 
the electron microscope back to the criterion for visibility in the 
light microscope. 


Contrast in the optical and electron microscopes. The term 
‘‘contrast’’ deserves some clarification. If a histological section 


is viewed in a light microscope, only a portion of the section will 
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be in sharp focus. Within this in-focus layer some tissue elements 
appear more or less stained and some appear not to be stained at 
all. The contrast within this layer refers to the different intensi- 
ties of light leaving the layer; that is, the principal cause of the 
observed contrast will be absorption within the in-focus layer. In 
other words, the ‘‘in-focus contrast’’ will not be very different if 
sections equal in thickness to the depth of focus are used or if 
sections much thicker are used. This may be checked by actual 
observation. Sections of tissue approximately equal in thickness 
to the depth of focus of a light microscope objective and stained 
with an organic dye, when viewed under oil immersion, show a con- 
trast which appears little different from the contrast seen in thicker 
sections. Actually, the thin sections appear clearer due to the 
elimination of superimposed out-of-focus images. 

As sections successively thinner than the depth of focus are ob- 
served, the contrast decreases, but only because there is less 
stained material, and even very thin sections of the order of 500 A 
show sufficient contrast for easy visibility. Sections of tissue 
stained with the (blue) basic dye gallocyanin chromalum and cut 
thin enough for electron microscopic observation have been ex- 
amined under oil immersion in the white light microscope. Baso- 
philic components appear blue, sometimes dark blue. 

The depth of focus of a white light microscope is approximately 
0.5u. Let us say that an electron microscope section is about 500 
A thick or about one-tenth the thickness of the depth of focus of a 
white light microscope. Unless indicated otherwise, it will be as- 
sumed that whenever optical sections are considered these sec- 
tions have a thickness of about half the depth of focus or about 
five times the thickness of an electron optical section. It is clear 
that such sections will have good contrast and cytological details 
will be easily visible. 

Let us consider a section cut for the optical microscope and 
stained with one of the common histological dyes. It will be as- 
sumed that monochromatic light is used for viewing and that Beer’s 
law is valid for light absorption by the dye. Then, for a concen- 
tration c of dye in an optical path length J, the light absorption 
will be governed by the equation 


] = [yer eto, (1) 
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where /, = intensity of light incident on that part of the section 
containing the dye of concentration c, 
] = intensity of light emerging from that part of the section 
containing the dye of concentration c, 
¢ = extinction coefficient. 
Let EF, be called the optical extinction, so that 


oe =€Cl). 


Let us consider a section cut for the electron microscope, the 
tissue having been stained in a similar fashion to that viewed op- 
tically. If a sufficiently small objective aperture is used, the 
scattering of electrons by the dye alone will be governed at least 
approximately by 


ey Poe ee (2) 


where /, = incident electron intensity, 
/ = intensity after penetrating a path length /., 
p = density of the dye, 
S = mass cross section. 

The value of § will vary with the accelerating voltage used. Its 
order of magnitude may be estimated from the literature. 

Lenz (1954) computed scattering cross sections for several ele- 
ments. He found good agreement between his computed angular 
distribution for small angle scattering and the experimental measure- 
ment of this scattering by Biberman et al. (1949). Lenz’s value of 
S was the order of 10°cm?/gm for accelerating potentials of inter- 
est to electron microscopy. Hall (1951, 1953, 1955), using a phe- 
nomenological theory of contrast production in an electron micro- 
scope with no objective aperture, reported S values of about 
3.5 x 10*cm?/gm for accelerating potentials of 65 kv. In this 
paper calculations using both values of S are presented. It will be 
seen that the qualitative conclusions are the same with either value. 

The choice of S as a constant independent of the chemical com- 
position is an approximation based upon two reports: (1) the report 
of Hall (1955) that S values were relatively constant despite wide 
variations in the atomic weight of scatterers used; (2) the curves 
of Lenz (loc. cit.) which indicate equal S values for chromium and 
gold and a very similar value for carbon. 

The quantity 7, will refer to the thickness of the section used in 
the electron microscope or the thickness of the tissue component 
under discussion, whichever is smaller. 
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Let EF. be called the electron extinction of the dye, so that 
E,=pSl,. 
If c is measured in moles per liter, then 
p = 107 °Me, (3) 
where M is the molecular weight of the dye. 


Comparing the two types of extinctions, we find the relation 


E St M 1 
cae eee Sig Se (4) 
E, €¢l, En Le, 
Assuming Lenz’s value for S, we find 
E Mt ec 
€_40?-— = —., (5) 
0 ely Cy 


If we consider a tissue component that extends throughout the 
thickness of the electron microscope section, then, under the above 
assumptions, we have 


L, = 5l,. (6) 


Values of « vary from about 10* to about 105 (Michaelis, 1947). 
As one example let us consider a dye of molecular weight about 
300, say, toluidine blue, viewed at a wave length at which e = 30,000. 
Then 


c 
= 0.20 -=. (Ta) 


E 
pee 
0 Co 


Using Hall’s value for the scattering cross section, we may write 


Be _ 0.07 <s (7b) 
at = 0. . 


0 Co 


According to Einarson (1951), staining by gallocyanin chrom alum 
occurs by the union of a lake-cation with negative valencies of 
nucleoprotein. However, regardless of the exact mode of binding, 
basophilism involves the addition of molecules having a molecular 
weight of about 400 to basophilic cell structures. N. Smith (per- 
sonal communication) has shown that if one accepts Einarson’s 
theory, « varies from about 10,000 to about 45,000. If we consider 
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a section stained with this dye and viewed so that « = 20,000, we , 
have 


E Cc 
—£=0.40 —. (8a) 
0 9 
Using Hall’s value, we may write 
E c 
ek ea Sp Oe, (8b) 
Ey Co 


As another example, the periodic acid-leucofuchsin staining of 
glycogen (Bondareff, loc. cit.) will be considered. 

Ultrathin sections of guinea pig liver, stained by Bondareff’s 
procedure, were examined under oil immersion in a light micro- 
scope with a Farrand interference filter No. 2091. This filter has 
a peak at 476 mu and a half band width of 17 mu. With this filter 
the fuchsin stained glycogen granules were easily visible. 

Stowell (1943) studied the absorption of light by a formaldehyde- 
fuchsin-sulfurous acid solution. From the data of Stowell the 
value of ¢ at 476 mu is approximately 30,000. 

Basic fuchsin is a mixture of pararosanalin and rosanilin having 
molecular weights of 324 and 337, respectively. 

Since the above observation was made on ultrathin sections, we 
may now assume 


l. = Los 
and we have 
E 830 ¢ c 
pig Abbe —=1.1-£, 
E 30,000 ¢ C ss) 


—* = 0.38 —£. 
- (9b) 


At first glance it might be thought that c, = c,, since both quan- 
tities refer to the concentration of dye. However, as will be shown 
in the next section, c, will be considerably greater than cy be- 


cause the resolution of the electron microscope greatly exceeds 
that of the light microscope. 


Resolution-density effect. In only a limited number of cases 
will a tissue component stain uniformly over dimensions equal to 
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the limit of resolution of the light microscope, say, 0.2u. In the 
light microscope all dyed structures will appear uniform over di- 
mensions of 0.2u or some dimension close to 0.24. However, this 
is the result of a smearing-out effect of the finite resolution fol- 
lowing the absorption of light by dyes distributed in submicroscopic 
patterns. This means that c,, the concentration effective in the 
electron microscope, will almost always be greater than cg, the 
concentration effective in the light microscope. 

The resolution-density effect is a familiar one in light micro- 
scopy. It will often happen that in a tissue section which shows 
little or no color to the naked eye, certain cytological structures 
will be intensely stained when viewed under oil immersion with 
the light microscope. 

If the stained tissue component extends throughout the thick- 
ness of the electron microscope section, equations (7) and (8) 
show that only a small submicroscopic concentration will make 
the electron extinction due to the dye alone comparable or supe- 
rior to the optical extinction. If the stained tissue component is 
thinner than the electron microscope section, equation (8) can no 
longer be used. In this latter case the order of magnitude of the 
electron extinction may be seen by taking specific examples. 

Let us suppose we are dealing with rat liver stained with gallo- 
cyanin chromalum. Finck (loc. cit.) found that the basophilia ap- 
peared largely, if not exclusively, in concentrations of the order 
of 200 A. Based on his work, estimates of the fractional volume 
occupied by such concentrates can be made. Estimates indicate 
that in a region that appears basophilic in the white light micro- 
scope only about one thirtieth of the region may contain the baso- 
philic granules. Hence, a conservative estimate would be 


Cc 
—£ = 10. 
Co 


For such a concentrate /, = 0.02 and, from equation (5), for 
gallocyanin chromalum, assuming that 1, = 0.25y, we find 


E 400 0.02 


a 82010" eee Wane] 6G. (10a) 
ES 20,000 0.25 
Using Hall’s value for S, we obtain 
E 
—£ = 0.56 (10b) 
Ey 
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Let us now consider two specimens, specimen one having been 
stained and embedded, while specimen two has been embedded 
without staining. Both specimens are sectioned identically and 
the extinctions of an electron beam traversing each section is ex- 
amined. The extinction of a structure in section one is due to the 
structure itself plus the dye plus any methacrylate which has in- 
filtrated the dyed structure. The extinction in section two is that 
due to the structure itself plus any infiltrated methacrylate. 

Let E, = extinction of the dry structure, 


E’, = extinction of the dye, 
E.1 = extinction of the methacrylate in the dyed structure, 
m2 = extinction of the methacrylate in the undyed structure, 


E', = optical extinction. 
It will be assumed that the extinctions are additive and hence 


Es =E. +E, +E 4; 
Ee = E, + E 2 . 
For a detectable contrast difference we must have 
E, . E, 2 q; (11) 


where g is a number that is small compared with unity. g is equal 
to the fractional intensity difference that is just detectable and is, 
for our purposes, primarily a psychological parameter. Chapanis, 
Garner, and Morgan (1949) report that under suitable conditions 
fractional intensity differences equal to 0.02 may be detected. We 
shall assume, somewhat arbitrarily, that g is equal to 0.05. Re- 
calling that we are neglecting contrast intensifications that may 
be achieved by using high contrast film it can be seen that g = 0.05 
is a conservative value. 
From equation (11) 


E, +E, ~ En, & gs 


OF Eg se +E a) 
Qqt+E on ee 
> m m 
0 ee fe 
Now ES =Sp i, rs 


where p,,, is the density of the methacrylate within the structure; 
and 


ey a 
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where p_ = density of methacrylate, 
¥, = fractional volume occupied by methacrylate within the 
unstained structure. 
Similarly 


where Pay =% 
m 


so that equation (12) becomes 


> q+Sp_,t. (vg - 2;) 
(E/E 4) 

Let us now consider the basophilic granules described by Finck. 
We have, assuming Lenz’s value for S, 


- 0.05 +10°x 1x 200 x 1078 (v, - 2) 
ue 1.6 
E , 2 0.03 + 0.18 (v, — v,). (14a) 


J (13) 


Assuming Hall’s value, we have 
E , 2 0.09 + 0.07 (v, — v,). (14b) 


The inequalities (14a) and (14b) are rather peculiar visibility 
criteria. They tell us what extinction must be produced in the 
light microscope in order that appreciable contrast differences may 
be achieved in the electron microscope. We must now see if we 
may reasonably expect the above inequalities to be fulfilled. 

We first note that v, — v, is certainly smaller than unity and, 
for the granules under consideration, perhaps quite a bit smaller 
than unity. 

We next recall what is meant by the term “‘a good optical dye.” 
This means, among other things, that tissues stain clearly and in- 
dubitably with the dye. In other words, the extinctions in stained 
tissue sections must be at least an order of magnitude above 


threshhold. In our example we therefore have 


E> 10 x 0.05 = 0.5, 


if we adopt the same threshhold value for light microscope obser- 
vations and if we choose a factor of ten to represent an order of 
magnitude difference. 

It is clear that either visibility criterion (14a) or (14b) will be 
satisfied and satisfied with a safety factor of about 5 or better. 
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It therefore follows that at least as far as this stain is concerned 
we are not dealing with a threshhold phenomenon but with a phe- 
nomenon that should be readily apparent. 

If we now turn to the leuchofuchsin staining of glycogen re- 
ported by Bondareff (loc. cit.) we have from equations (13), (9a), 
and (9b) 


_ 0.05 + 0.5 (v7, - 24) 


One ’ 
tbe fe, 

0.05 + 0.18 (v, - v,) 

or 2 : 
: 0.38 ¢./c, 


It is clear from the previous discussion that, even with values 
of c,/c, that are not too high we may expect the above criteria to 
be fulfilled. 

We now attempt to generalize from the above cases. We note 
that our qualitative conclusions are not too sensitive to a reason- 
able choice of S nor to variations in 7, that are not too large. We 
next note that the factor E,/E, will, apart from the resolution 
density effect, not vary too significantly from one dye to another 
since dyes commonly used have molecular’ weights of, roughly, 
250 to 1000. We may therefore expect a wide variety of classical 
organic dyes to be useful for electron microscopy. 

It may happen that upon comparing a stained section with an un- 
stained control the stained section will show less contrast. This 
will occur if the tissue elements that are stained have a lower 
density than their surroundings and if the staining is not strong 
enough to reverse the contrast. Such a result would, of course, be 
a special case of the more general conclusions given above. 

The above considerations may be generalized still further. The 
comparison of electron extinctions to optical extinctions treated 
above was simply a stratagem to demonstrate the usability of or- 
ganic dyes for electron microscopy. Actually, of course, electron 
extinction depends in no way upon color properties and hence we 
may expect a wide variety of organic substances to have useful- 
ness as electron stains. Two conditions appear necessary: (1) 
The electron dye must combine with certain tissue elements in 
high concentrations. (2) The molecular density of the dye must be 
higher than methacrylate, the higher the better. 

It may be expected that cytologists will eventually develop a 
variety of electron microscope dyes of known cytological proper- 
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ties which satisfy the second requirement above by incorporating 
heavy atoms into organic molecules. However it may be seen that 
a wide variety of stains is already at hand which can be applied 
usefully to submicroscopic morphology and histochemistry. 


Effects of the electron beam. As demonstrated by Morgan, 
Moore, and Rose (1956), changes in contrast may occur in tissue 
sections due to interaction with the electron beam. These changes 
can occur at beam levels below those ordinarily used for observa- 
tion. Consequently, the contrast in the image that we ordinarily 
observe will not necessarily reflect the contrast that existed be- 
fore beam damage. The beam may increase or decrease the con- 
trast, and changes apparently depend on a variety of factors such 
as prior tissue treatment and quality of sectioning. How electron 
beam damage will modify the contrast in the cases considered 
here is as yet unknown. This effect may demand that quantitative 
verification of the contrast considerations developed be done only 
at beam levels below those currently employed in electron micro- 
scopes. However, in view of the evidence that we are not dealing 
with threshhold phenomena we may proceed, subject to later verl- 
fication, under the assumption that qualitative conclusions may 
still be drawn by staining with organic dyes. 


Miscellaneous considerations. Sections may show appreciable 
contrast in the electron microscope while showing no contrast at 
all in the light microscope. Unstained specimens and metal- 
stained specimens are examples. An analysis of such cases in- 
volves not a comparison between optical and electron contrast, 
but an absolute calculation of contrast in the electron microscope 
alone. In general, such a calculation appears to be impossible. 
There is no way of knowing, @ priori, the distribution of any spe- 
cific electron microscope dye, nor is there any @ priori way of 
knowing the distribution of protein or nucleoprotein in unstained 
tissue. Consequently, there is no way of predicting the contrast 
in any specific case. However, it is possible, and may be useful, 
to make such estimates @ posteriori for the purpose of showing 
that given results are reasonable. 


Unstained tissue. For example, Gersh et al. (loc. cit.) found 
that submicroscopic vacuoles could be obverved in unstained sec- 
tions of guinea pig liver. Is it reasonable that unstained material 
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should show sufficient variation in density so that such vacuoles 
might be observed? 
As in equation (2) 
I =1,e7?8'e. 
Let /, be the intensity after passing through a wall, /, the inten- 


sity after passing through a vacuole. Let p, and p, refer to the 
respective densities, 


ur — gle (P1~ Pa). 
I, 
Let S = 105 cm?/gm, 
1, = 500 A =5 x 10-% cm. 
Then 


I, 
P43 — Pg =-2 a 
2 
If we again somewhat arbitrarily choose g = 0.05 and hence /,//,< 
0.95 for a criterion of visibility, then we obtain 


P1—P, 2 9-1. 
The method of tissue preparation described in Gersh et al. (loc. 
cit.) should extract the tissue water and some lipids, replacing 


these by methacrylate. Since methacrylate has a density of about 
1 gm/cm®, we may approximate Po by 


p. = 1 gm/cm’®. 


Actually p, will be a little larger, since the vacuoles contain some 
protein. However, if we ignore this we have 


p, 21.10 gm/cm?. 


Finck found nucleic acids of liver cells mainly distributed in 
the walls of submicroscopic vacuoles. Nucleic acid has a density 
of about 1.6 gm/cm*®. The proportion of nucleic acid, protein, and 
methacrylate in the walls is unknown, but in view of the high con- 
centration of nucleic acid in the walls it is certainly not unreason- 
able to expect that the relationship 


p, 2 1.10 gm/cm? 


will be satisfied, and hence it is not unreasonable that, though un- 
stained, submicroscopic vacuoles of liver cells should be visible. 
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It should be mentioned, however, that submicroscopic vacuoles 
may be seen following ribonuclease digestion even though the con- 
trast between walls and vacuoles is greatly reduced. This is in- 
terpreted to mean that the protein component of liver cells is more 
concentrated in the walls than it is in the vacuoles. 


Enzyme localization. Some preliminary work by L. Nelson in- 
dicates that it may be possible to localize the sites of certain 
enzymic activities at an electron microscope level. The possibil- 
ity of localizing enzymes on a submicroscopic level is of obvious 
importance. We shall not now consider the difficult questions of 
specificity and diffusion artifacts, but simply indicate that con- 
siderations of contrast should not serve as a limiting factor in 
such studies. 

Doyle, Omoto, and Doyle (1951) studied phosphatase activity in 
kidney cortex and reported as much as 0.1u mol of PO,= deposited 
in a tissue section 10 mm? x 5u. This is equivalent to a density 
of Ca,(PO,), of about 0.1 gm/cm®. From calculations above it is 
clear that observable contrasts will be obtained with the electron 
microscope even if much smaller amounts of insoluble calcium 
salts were deposited, especially since the deposition will un- 
doubtedly show segregation in submicroscopic structures. 


SUMMARY. The extinctions produced by certain organic dyes 
in the light microscope and the electron microscope are compared. 
It is shown that the two have the same order of magnitude. The in- 
fluence of the methacrylate embedding medium is not sufficiently 
great so as to hamper visibility in the electron microscope. It 
is pointed out that a wide variety of organic stains should be use- 
ful for electron microscopic studies. 


The author wishes to express his great appreciation to Dr. Isi- 
dore Gersh for many helpful discussions throughout the course of 
this work. This work was supported in part by grants from the 
Wallace C. and Clara A. Abbott Memorial Fund of The University 
of Chicago and from the Commonwealth Fund of New York. 
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In electrocardiography the electrical potentials due to the heart ac- 
tions can be analyzed by assuming the human body to be a conductor of 
homogeneous medium and the heart to be a combination of singularities 
within it. For a spherical conductor the ‘ ‘interior sphere theorem’? of 
G. Ludford, J. Martinek, and G. Yeh (Proc. Cambridge Phit. SOC, y ol, 
389-93, 1955) renders potential expressions due to any singularity. For 
a conductor of prolate spheroidal shape the potential expressions due to 
a source-sink pair and a general dipole have been given by J. R. Wait 
(Jour. App. Physics, 24, 496-97, 1953) and the authors (paper at the 
Conference on the Electrophysiology of the Heart, Feb. 16-17, 1956, in 
New York, to appear in the Ann. N. Y. Acad, Sciences) respectively. (A 
theorem which applies to any singularity inside a prolate or oblate sphe- 
roid will be published by the authors soon.) This paper presents numeri- 
cal and graphical results of potentials on the surfaces of a prolate sphe- 
roid and a sphere due to source-sink pairs and dipoles of several loca- 
tions and directions and compares the various representations. A dis- 
cussion on the judicious choice of heart models concludes the paper. 


Introduction. Some new contributions to the mathematical 
theory of electrocardiography have been reviewed recently (Martinek 
and Yeh, 1956). In particular, several theorems deal with the elec- 
trical potential due to the heart actions by assuming the human 
trunk to be a spherical or spheroidal conductor of homogeneous 
medium and the heart to be a source-sink pair or a general dipole. 
It is the purpose of this paper to demonstrate the evaluation of the 
surface potentials by means of these theorems and to compare the 
results due to various representations of the heart and the torso. 


*This investigation was supported by The National Heart Institute 
under a research grant H-2263. 
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A general dipole in a prolate spheroid. Consider a prolate sphe- 
roid of semi-major axis a and semi-minor axis }=ay1- e” where 
0 < e<1is the eccentricity. The equation of its surface in Carte- 
sian coordinates is given by 


ew? + y? 


Se, ehet ea 5 (1) 


A dipole of strength » and direction along the unit vector B 
B=B,i+B,j+B,k=Bz,a,+B,a, +B a, (2) 
is located at the point (z,y,,2,)- In equation (2) the prolate 
spheroidal coordinates (€, 7, ~) have been used: 
ee 
# = ae[(1-7”)(€? — 1)} cos 9, 
ae 
y = ael(1- 9*)(€? - 1}? sin 9, (3) 
2= ae En ’ 


with 1<&€0, —-1€7y <1, and 0S @€<2z. The components 


Sere , and B., can be expressed in terms of the components B.; 
By, seh B, as follows: 


p, = Be 22 , B, ay B, oe 
= he 0€ he 0€ he 0€ 
Boa B, a i og 
pee e ae Se (4) 
h, on h,, an h,, on 
B im dx By ay OF 
Glin SIS Ma oe oe ’ 
hg OP hy ap hy OP 
where 
€2_ 7? 
he = ae pins 


Tre (5) 
Ree eee 
1-y 


hg = aeV(E?— i)(i- 7). 
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In G. Yeh and J. Martinek (1956) the potential on the surface of 
the spheroid € = €, is found to be 


&(E.7,?)= > D. (47, cos m(P - Pq) + AZ, sin m(P - Po)I + 


n=1 m0 


m al (é,) m m 
es (€,) rT P™ (é,) ee (e0| sea (n), (6) 
where 
2 
Ar =i (20+ Ie” (-1)™ foam . 
ae (n + m)! 
Be, , B,, , 
5 Pat (9) Pr (£0) + > Par (0) Par (£0) » (7) 
€o ul) 
~ .# Boo a[(m— ml? m 
A’ = ac (2n + 1) me, (-1) C= gal CES) ek CY 


In Eqs. (6) and (7) the following notations are used: 

(1)«, =1form=0 ande, =2 for m #0. 

(2) P™(é) and Q® (€) are associated Legendre functions for 
arguments greater than one as defined by the Mathematical Tables 
Project (1945). The expressions (7) differ from those in Yeh and 
Martinek (Joc. cit.) by a factor of 7” where these functions are de- 
fined according to P. Morse and H. Feshbach (1953) which are 
somewhat different from those defined by the Mathematical Tables 
Project. 

(3) P™(n) = P™(cos 6) with 0 <€6<90° is the associated Leg- 
endre function for argument less than one as defined by the Mathe- 
matical Tables Project (Joc. cit.). Note that for negative arguments 
P™(—n) = (-1)"*" P™(n). 

(4) The derivatives are 


Z d 
ee | 


: d 
Qn (61) = ee Qn a 


dé 2 Bp (8) 


P™ (cos 0) = 
Ox 


99 


= d 
Bee EB P™ (cos o| 
sin 4, G0 = 6, 


Benge SR STN F cos 0 
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Numerical example I. A numerical example has been carried out 
for a prolate spheroid 6 = 2 aor e= V5/3 = 0.745 with a dipole of 


strength pat 2% = 2%, = 2/3, Yo = Yo, = 2/6, and Zy = 29,= 4/4 


pointing towards the point 2, = @/3, Yo. = — a/6, and 2), = Ta/12; 
that is, 
Bae cf? erat ae 
2 a 
[(@ gq ~ 1)? + (Yoo ~ Yor) + (202 ~ 201) 4 
Yoo — Y 
B, 3 2 oS eee — = -0.707, (9) 


kf Zo)" pe Joi)” + (2e5 = a el st 


2 ee 
hed tet ane Eas SOR a ee ee Oe 


ort 
leper 1)" + (Yoo - Yor)” + (295 — 204 r 


From equations (3) we have 


2 2 2 2 
fo 1 Pao tos Te ot: 4 
wo) itace Gee eee Whee 
r 2 2,2 


02, ae ‘ 
ya tS 1.2742 
ft @o1 + Yor + 201 ol 42o,| | _ ; (10) 
a*e? ae? 0.0883 
Yor _ rf 
Po, = arc tan —— = are tan >. 
Zo, 


Hence we obtain €,, = 1.1288, ,, =0.2972=cos 72°43’ and 


Py = 26°34’. . 
From equations (3), (4), (5), and (9) we obtain 


Be By Boo: 
ae ——~ = — 0.1019, ae --o = 0.6889. Ge = 1.2647. (11) 
h h 
fo1 No1 Po1 


On the surface €=&, = , 


= 1.34, n = 2/aeé, = 2/a and ~ = 


arc tan .. The computations have been performed for the surface 


potential at a number of points. The numerical results by taking 
m= 1, 2, and 3 in the series equation (6) are shown in Tables 1 to 
5, and the graphical results are shown as the curves I in Figures 
1 to 5. The curves indicate positive potential values on the outer 


. ‘a teatiediat | 


Wn a nal a he tia A 
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normals at the surface points and negative potential values on the 
inner normals. 


TABLE 1. Surface potentials at 7 =0 orz=0. 


I 


=] -0.647 |-0.535 | -0.248 - 0.248 
- 3/4 - 1.821 |-0.916 | -0.136 
-1/2 -2.254 |- 1.590 0.00 1 
-1/4 -1.779 |-1.742 0.290 
- 1/8 - 1.030 |-1.576 0.489 
-1/16 -0.643 |- 1.429 0.601 
0 -0.217 |-1.240 0.718 
1/16 0.245 |-1.011 0.841 
1/8 0.735 |-0.744 0.968 
1/4 1.774 |-0.102 1.228 
WS, 3.858 1.496 1715 
3/4 5.340 3.152 1.999 
1 8.614 2.907 1.333 
eS 
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TABLE 3. Surface potentials at 9 =0° and M = 180°. 


180° 


TABLE 4. ; ~ = 30° 
TABLE ‘| Surface potentials at ; Sagal 


Imagine the prolate spheroid to be a human trunk with the z-axis 


pointing to the front, the y-axis pointing to the left, and the z-axis 


pointing to the head. Then the surface potentials computed are 
those due to the heart actions represented by a dipole of strength 
w located at a front-left-upper point, pointing 45° towards upper- 
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right. The body shapes and singularity locations are shown in 
Figure 6. ; 


A source-sink pair in a prolate spheroid. If the heart actions are 
represented by a source of strength x at the point (x, y¥,, 24) or 
(€,, 4, 94) together with a sink of strength — x at the point (a_, y_, 2_) 
or (€_,_,P_) instead of a general dipole, ‘the surface potential is 
given by J. R. Wait (1958) as follows: 


®(E,mP)= DY, (Ang COS ™(P — P+) — Bag CO8 m(P - G_)]+ 


n=1 m=O 


, (12) 
Qn (é,) 
Cee Pe ee.) 
| te oe lea 
where 
x (n- m)!\? 
A. = —(2 1 —1)™ |————_|_ P™ oe 
ks a n+ le, (-1) een Gs) Lh Sa) 5 
(13) 
x erty eo 
B= —(2n+ De—t{—3) BEG PE CEs) 2 
ae (n + m)! 
Two numerical examples have been carried out for this case. 
Numerical example II. Take aw = 2, = Al 3, yo. 2°¥5, = a7 SG; 


Z_ = By, = @/4, and a = &q_ = A/3, Y4 = Yoo = — a/6, 24 = 29, = 10/12; 
that is, the sink is located at the same point as the dipole in Ex- 
ample I and the direction from the sink to the source is the same 
as the direction of the dipole in Example I, the source being located 
at a front-right-upper point. Then from equations (11) we have 
é_ = 1.1288, n_ = cos 72°43’, and 9_ = 26°34’; and similarly we 
have &, = 1.1916, 7, =cos 48°57’ and 9, = — 26°34’. The nu- 
merical results are summarized in Tables 1 to 5, and the graphical 
results are shown as the curves II in Figures 1 to 5. 


Numerical example Ila. Take 2 = 2, = a/3, Y~ = Yo, = 2/6, 
Z.= 2, = @/4, and a, = Zo, = a/ 38, Ys = Vox =% 24 = 203 = 5a/ 12; 
that is, the sink is located at the same point as the dipole in Ex- 
ample I and the direction from the sink to the source is the same 
as the direction of the dipole in Example I, the sink being located 
at a front-center-upper point. Then from equation (11) we have 
é_ = 1.1288, n_ = cos 72°43’, and 9_ = 26 °34’; and similarly we 
have €, = 1.1180, 7, = cos 60°, and 9, =0°. The numerical re- 
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FIGURE 1. Surface potentials at n = 0 or 2 =0. 


sults are summarized in Tables 1 to 5, and the graphical results 
are shown as the curves IIa in Figures 1 to 5. 


A general dipole in a sphere. Consider a sphere of radius 6 = 
a with its center at the origin. This is an inscribed sphere in 
the prolate spheroid in Examples I, I, and Ia. 

In G. Ludford et al. (1955) an ‘‘interior sphere theorem’ is 
established which states that if ®)(r) = ®)(r,0,p)(7,6,9 are 
spherical polar coordinates) is the potential of any singularity dis- 
tribution of total mass zero in r < 6 in an infinite domain, then 


2 co 2 
®(r) = ©(1, 0, 0) = ©, (r) + : ®, ().2 f A®, (AS) ds (14) 
b 


r br r 


is the potential due to the same singularity distribution in the 
sphere. 
For the specific singularity of a general dipole M at the point 


(with w < 6) the potential expression has been derived in G. Yeh 
et al. (1955). 
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b? 
8 | | ere 
(erence 8 a Se ( oe ) 


Sites lr—o[? 52 3 
oii a 
of 
42 (15) 
r2- —(@-t 
M,-ft w-f me ) M, 

-r\? b? be’ 
@ w 


where M, = (M-@/o*)@ is parallel to w and M, =M-M;, is per 
pendicular to w (boldface denotes a vector) and r=2i+ yj t+ 2k. 

Two numerical examples have been carried out for the surface 
potential on r= 6 and z= 0. 


Numerical example 1, Taking the same dipole as in Example I, 
we have M = nB = 0.707p:(-j +k) and @ = a(i/3 + j/6 + k/4). Then 
M, =1(0.0982i + 0.0491j + 0.0737k), M, = n(- 0.09821 - 0.756] + 
0.633k), M =p, and w =0.4488a. The potential expression, equa 


FIGURE la Surface potentials at 7 = 0 or 2 =,0. 
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y 


FIGURE 2. Surface potentials at 7 = + or z= "eh 


tion (15), on r= 5 and 2 =0 with z= bcos 9 and y= 3 sin @ be- 
comes 
a? e* 0(8) — 0.8838 (sin + 0.1259) 
Se 
lL [(cos p — 0.5)? + (sin p — 0.25)? + 0.1406]? 


0.8132 (cos m + 4.099 sin p — 1.015) 
[(cos p — 1.111)? + (sin p — 0.5556)? + 0.6944]? 
0.1824(cos p — 7.699 sin 9) 
1 - 0.1889(2 cos 9 + sin 9)? 
1 - 0.5556(2 cos 9 + sin 9) 


: (16) 


0.3714 (2 cos 9 + sin ) — 0.2456. 


The numerical results are summarized in Table 1, and the graphical 
results are shown as the curve I, in Figure 1a. 


————=— Fhe NY 


Se AS 
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Numerical Example I,,. Take at the same point w = a(i/3 + j/6 + 
k/4) a horizontal dipole M = — ,j pointing from left to the right. 
Then M, = — 2(0.27781 + 0.1889] + 0.2083k), M, = n(+ 0.27781 = 
0.8611j +.0.2088k), M=p and w»=0.4488a The potential ex- 
pression, equation (15) on r= 6 and zg = 0 becomes 


a’ e? 0 (8) ~1.25(sin 9 ~ 0.2257) 
= se 
p [(cos @ — 0.5)? + (sin » ~ 0.25)? + 0.1406]? 
2.276 (cos o - 1.300 sin 9 - 1.014) 
4+ 


[(cos o ~ 1.111)? + (sin 9 — 0.5556)? + 0.6944]? 
0.5158(cos » — 3.100 sin 9) 
1 — 0.1389(2 cos 9+ sin 9)? _ 


(17) 


1 — 0.5556(2 cos 9 + sin 9) 
+ 
Hass p — 1.111)? + (sin p - 0.5556)? + 0.694432 
0.38714(2 cos 9 + sin a) — 0.6944. 


The numerical results are summarized in Table 1, and the graphical 
results are shown as the curve I, in Figure 1a. 


A=0° plane | 
$ =180° 7 


FIGURE 3. Surface potentials at p = 0 ° and ~ = 180°. 
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ol 


FIGURE 4. Surface potentials at P = + 30°. (Top) 
FIGURE 5. Surface potentials at 9 = — 30°. (Bottom) 


A source-sink pair in a sphere. Another special case of equa- 
tion (14) is the potential due to a source of strength x at the point 
w, together with a sink of strength —x at the point w_ which has 
been derived in G. Yeh et al. (loc. cit.) as follows: 


% % x 
© (9) = ——— = ee ee 
eel eee eed @ - e.| 
oS OPS fees Wiel ie 
(18) 
b? 
b? - w_ r+o.|r- 7 o_| 
x - 
7 108. p2 
b —@y'Pt+ oy, Sale 
- 


Three numerical examples have been carried out for the surface 
potential on r= 6 and z = 0. 
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Numerical Example II,. Taking the same source-sink pair as in 
Example II we have 


@— = Fb + Yo.) + 29, = a(i/3 + j/6 + k/4), 
and ; 
Wy = Boi + Yoo) + Zo, k = ali/3 - j/6 + TR/12). 


Then w_ = 0.4488 a and w, = 0.6923 4° The potential expression, 
equation (18) on r= b and 2=0 with x=) cos and y= 6 sin 
becomes 


1 
os ©(b) = -1.118[(cos p — 0.5)? + (sin p — 0.25)? + 0.1406] ? + 
1.118[(cos p — 0.5)? + (sin p + 0.25)? + 0.7656)" ? - 


1 
2.491[(cos @ — 1.103)? + (sin p — 0.5516)? + 0.6848] 7 + 


Sa, £19) 
1.615 [(cos 9 — 0.4636)? + (sin 9 + 0.2318)? + 0.1209] * + 


1.118 log, {1 - 0.5 cos p — 0.25 sin p + 0.6732.[(cos p — 1.103)* + 
(sin p — 0.5516)? + 0.6848]? }/ {1 ~ 0.5 cos p + 0.25 sin 9 + 
1.038 [(cos o — 0.4636)? + (sin @ + 0.2318)? + 0.1209} }. 


The numerical results are summarized in Table 1, and the graphical 
results are shown as the curve II, in Figure la 


Numerical Example II,,. Taking the same source-sink pair as 
in Example II, we have 


@_ = 2,1 + Voy) + 20, Rk = a(i/3 + j/6 + k/4), 
and Oy, out ti Vog ht Zogk = a(i/3 + 5k/12). 


Then w_ = 0.4488 a and wo, = 0.5336a. The potential expression, 
ee Se (18) on r= 6 and z = 0 becomes: 


~ (5) =—1.118[(cos » — 0.5)? + (sin p - 0. ee + 0.1406] z ip 


1.118[(cos p ~ 0.5)? + sin? p +0. 3906]? 


1 


2.491[(cos p — 1.108)? + (sin p — 0.5516)? + 0.6848] * 
2.095 [(cos p — 0.7804)? + sin? p + 0.3426)” 7 4 ay 
1.118 log, {1 - 0.5 cos g - 0.25 sin p + 0.6732 [(cos  - 1.103)? + 
(sin p — 0.5516)? + 0.6848]? }/ {1 - 0.5 cos @ + 


0.8004 [(cos p — 0.7804)? + sin? 9 + 0.3426]7} . 
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The numerical results are summarized in Table 1 and the graphical 


results are shown as the curve II, in Figure la. 


Numerical Example Il,,. Take the sink at the same point as in 
Examples II, and II,,. w_ = a(i/3 +j/6 + k/4) and the source at 
the point w, = a(i/3 - j/6 +k/4); that is, the direction from the 
sink to the source is the same as the direction of the dipole in Ex- 
ample I, (horizontal). The potential expression, equation (18), on 
r= 6 and zg =0 becomes 


eds 
% (bd) = — 1.118[(cos o — 0.5)? + (sin-p — 0.25)? + 0.1406] ? + 
m6 


ah 
1.118[(cos p ~ 0.5)? + (sin p + 0.25)? + 0.1406] ? - 


set 
2.491[(cos m — 1.103)? + (sin p — 0.5516)? + 0.6848] ? * (93) 
1 


2.491[(cos » — 1.103)? + (sin @ + 0.5516)? + 0.6848] 7 + 
1.118 log, {1 - 0.5 cos p — 0.25 sin @ + 0.6732 [(cos p - 1.103)? + 
(sin p — 0.5516)? + 0.6848]? t/{1 — 0.5 cos o + 0.25 sin 9 + 
0.6732 [(cos » — 1.103)? + (sin @ + 0.5516)? + 0.6848]7t. 


The numerical results are summarized in Table 1, and the graphical 
results are shown as the curve II, | in Figure la. 


CONCLUSIONS, Tables 1 to 5 and Figures 1 to 5 show the com- 
parison of surface potentials due to various representations of the 
heart and the torso shape. It is seen that a great variety of surface 
potentials can result by assuming the heart to be either a general 
dipole or a source-sink pair. The position, the direction, and the 
distance between singularities all contribute to the difference. 
Even if the location of a sink and the direction from a sink to a 
source are the same, the different distances between them give 
quite different surface potentials, and a dipole cannot be approxi- 
mated by a source-sink pair unless placed very closely apart. In 
addition to these two models, ‘equation (14) allows computation of 
potentials on a spherical surface due to any other singularity dis- 
tributions such as a vortex ring, a quadrupole, an octopole, etc., 
and the line, surface, and space distributions of such singularities 
and fundamental singularities. A similar formula for potentials on 
a spheroidal surface due to other singularities will be published 
by the authors in the near future. The judicious choice of heart 
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Singularity Locations 
ExamplesI&Io: dipole at P; pointing to P2 
ExamplesII & Io'sink at P; and source at P2 
Examples Il, & Igo'sink at P, and source at Ps 
y Example I}, ‘dipole at P, pointing to Py 

Example II) 'sink at P, and source at Py 


P, is at (2,2,2 ) 


; ie 
ps at (9,-2.9) 


P; is at (9,0,29) 
p, is at (9,-2,4) 


(a)Top View 


z y— the spheroid 72 


/ 


iz the inscribed 


eee Oe 


right left front back 
(b) Front View (c) Left Side View 


FIGURE 6. Body shapes and singularity locations. 


models depends on more numerical and graphical results due to 
various singularities as demonstrated in this paper and compari- 
sons with measured potentials due to healthy and sick hearts. 
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In the first part of this paper we have assembled some properties of 
the quantities RR”, where Rf, denotes the number of distributions of n 
different objects into m indifferent parcels, with no empty parcels al- 
lowed. We then discuss the following problem (N. Rashevsky, 1954, 
1955 a,b, 1956): to find the total number, G,, of graphs that can be ob- 
tained from the biotopological transformation 4 DX) for a given value of 
the parameter n. This is related to the distribution of n indifferent ob- 
jects into m different boxes. A formula for G, is given which, however, 
is not very convenient for practical computations because it involves a 
summation over certain ‘‘admissible partitions’? of the number mn (m is 
a second parameter of the transformation). Some theorems are derived; 
with their help we can simplify the calculation of G, to a small extent. 
The numbers G, are calculated for n <9 and estimated for n= 10. Itis 
found that G7 ~ 5.4 x 104, Gg = 8.3 x 10°, Gg = 1.4% 107, and Gig ~ 3x 10°. 
These values of n are those which might be used in connection with 
N. Rashevsky’s work (cf. Rashevsky, 1956). 


N. Rashevsky (1955a) has introduced certain quantities R* in 
trying to determine the total number of transformed graphs resulting 
from his biotopological transformation (T X). In W. A. Whitworth’s 
or P. A. MacMahon’s terminology #? is the number of ways of dis - 
tributing n different objects into m indifferent parcels with no empty 
parcels allowed. (Whitworth, 1951, Chap. Ill, also available in 
1925 edition, N. Y.: G. E. Stechart & Co.; MacMahon 1915, Chap. 
II). We shall derive some properties of these numbers RP”, mostly 
in the form of recurrence equations. Practically all of these are 
found scattered in the literature, but often in a form which does not 
show at once the connection to the combinatorial problem in which 
we are interested. 
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Calling ae the number of distributions of our n different objects 
into m different parcels, we have 


eee ey (1) 
and 
Rn = > '(-1)" (7) (m~ oy (2) 
v=0 
Therefore: 
n = v (m - v)” 
RP = Lu Cee are aera (3) 


which is Rashevsky’s (1955a) equation (1). 

We shall put R® =0 if m=0 but n> 0, and RF” =1lifn=m=0. 
The right-hand side of (3) is zero for m>n, cf. equation (17) on 
page 249 of KE. Netto’s book (1927; also available in 1901 edition, 
Leipzig: B. G. Teubner). 

Equation (2) is usually derived by the method of ‘‘Ein-und 
Ausschaltung’? (EA-method). In fact, (m—v)" is the number of 
ways of distributing the n different objects into m different parcels 
if v specified of the parcels are left empty, and the remaining par- 
cels may or may not be empty. This then gives equation (2) as the 
result of a general principle. For details see: Th. Skolem, 1917, 
pp. 4 ff., or 1927, $$ 179-180, pp. 281-83; Whitworth, Joc. cit., 
Chap. HI, propositions XIV and XXII; cf. also: G. Pélya and G. 
Szeg6, 1954, Volume II, p. 119 (exercise 21), p. 326 and references 
given there. 

If we write the 


‘ 


‘generating function’’ 


G_ (a) =(e* -1)™ (4) 


in the form 


G,, (@) = “% (”) (=1) gees 


v=0 


and then calculate the derivatives of G_, (v) for 2 = 0, we see at 


once that 
1 142° G {2 
‘ies de x=0 
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Poa a 
Thus Fk” is a1 times the coefficient of 2” in the expansion of 


(e* — 1)™ (Whitworth’s proposition XXIII, loc. cit., Chap. III): 


oo 


m! 
(e205 =e — Fae. (6) 


n=0 


Rashevsky’s (1955a) equation (12), namely, 


ab n—m 
R" 2 a n n—A 
=1 


follows immediately from (5) if we use Leibnitz’s rule for the de- 
rivatives of a product. 
We have 


R™=1 (8) 
and 


Re 1 (9) 


n 


for all n. Equation (8) is obtained directly from (3), and (9) from 
the second of Netto’s equations (17) (loc. cit., p. 249). 
Furthermore, 


R® =mR™*+R™'. (10) 


This recurrence equation, which provides a simple method for 
the calculation of the numbers F”, can easily be verified by in- 
serting the values of R”, R®~', and R”~' [obtained from (3)] into 
(10). 

Equations (7) or (10), together with the ‘‘initial values”’ oy ze 
completely equivalent to (3) in defining the numbers fe tor 1s €n. 

As an immediate consequence of (10) we find that® 


nmm 
=): fae ieee. (11) 

v=0 
A table of &* is found on page 169 of Netto’s book (loc. cit.) for 
the first eight values of n. (He writes (ae: instead of R” , with his 
9 corresponding to ourn.) Netto also gives a combinatorial interpre- 


*For a direct proof of equation (11) see appendix 1. 
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tation of our equation (10) which is the same as his equation (1) 
on p. 169. 
We shall now prove that 


n! 
ao , = ii2 
2D fee tere, 


(Part.) 


where the non-negative integers a,, @,,---,@, specify a partition 
of n into m positive parts, i.e., 


Ct Pee wie ee Ps). (13) 
or 


m=G,:1+a,:2+...-+ A, °K with (14) 


M=A,+G,t+.-.+@,3 4,70, k=n-mrtI, 


and the sum in (12) has one term for each such partition. 

By a composition of the number n into m parts we mean a parti- 
tion of that number in which the order of occurrence of the parts is 
essential (MacMahon, loc. cit., page 150). We shall usually denote 
a composition by a sum of greek letters, e.g., 


Tas hy tel gt caret A ng (15) 


m 


it being now understood that each permutation of the parts in (15) 
must be counted as a different composition. Since 


is the number of compositions that can be formed from one given 
partition like (14), we may write equation (12) in the form 


1 n! 
| ein a 
a ist ea) D fo Ag! ws. si} 6) 


(comp.) 


where the sum now is over all compositions (15) with positive \’s. 
It is convenient to introduce the quantities r”, related to Rr by 
the equations 


(17) 
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It follows then from (7) that 


= v} 

r= Lae ae 18 
rape Ce las ae 
uL=0 

and from (8) that 
¥ 1 
ri = ' for all v- (19) 
vt+l 


Finally, equation (16) becomes 


v! 
“aD - +I)! (0, +1)... (, + 0! - 


(Comp.) 


the sum now being taken over all compositions of the form 


V=O,+Ggt--+t+ Fm, with o, 20. (21) 

After these rearrangements it is very easy to show by induction 
on m that the right-hand side of (20) satisfies the recurrence equa- 
tion (18). Moreover, it is equal to (v + 1)~1 for m=1, and this 
completes the proof which we have taken from N. E. Norlund (1924, 
pp. 139-40). 

A different proof of equation (12) consists in expanding the ex- 
ponential function of G, (a), 


a xn? ot m 
G (2) = wet ayers tre ; 


finding the coefficient of 2” by means of the multinomial theorem, 
and comparing it with the same coefficient in equation (6). See: 
E. Netto, loc. cit., page 185. 

The quantities r¥ of equation (17) are related to the so-called 
Bernouilli polynomials of negative order—more precisely, the 


Bemouilli polynomials with equal spans, @, = @_ =+++ = Om = i 
These polynomials, BO™ (2), may be defined by the equation 
eS | ot eae eee 5 
Kr (22) 


v=0 
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(Nédrlund, loc. cit., equation (75), p. 142, also p. 145). Putting 
2 = 0, and 


Shegana (1) gi priatae 


(22) becomes: 


foe) tn 
f= A)" = Ree Sey 5 ig 23 
G ) oe (n—omy ft 2™ So 
Comparison with (6) shows that 
R = (7) Bo™ (mn). (24) 


This equation allows us to derive some more relations between 
the numbers #”. For example, Nérlund’s equation (50) (loc. cit., 
p. 189) gives, fora = y=0,n ~+~-q <0, p > -p <0: 


Vv 

Dera + qty _ aig al! fa qts ppty-s 

ile sa eda FE, 
s=0 


or 
nies, — oe 
(7) R” = FES aL er (25) 
A=0 


in which q is any fixed integer satisfying 0 <q <m. 
This equation may also be written in the form 


rine X pn-v 
Ra)? (x) #e Be -q? 


and this can easily be verified directly, either by using the defini- 
tion of R, equation (2), or also with the help of equation (5). 
For ¢ = 1 or g = m — 1 equation (25) becomes the same as (7). 
Again, if in the first equation on page 146 of Nérlund’s book 
(Joc. cit.) we replace n +1 by -(m-1) <0 and n by -m, the re- 
sult, expressed in terms of the R’s, will be 


m Vv 
Rm-ity_ m+v m+v— 
mt (m+n) 2 8 ) 8. Rs - (26) 
s=0 
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where B. are the ordinary Bernouilli numbers. Now, in equation 
(26) write n instead of m+ v and replace the left-hand side by 
R® —~mk®-' [from equation (10)]. Remembering that By = 1, 
Bee = 53" and all other B, with odd index s are equal to zero, we 
obtain after some simple rearrangements 


(n — m) - Oe et — n n-s 
ren = Rts y (7) BR (27) 
s=2 


where the prime indicates that the summation may be taken over 
even values of s only.* In this recurrence equation the index m is 
the same throughout and only the upper index varies. 

For the sake of completeness we should also draw the reader’s 
attention to exercises 765-769 (p. 272) and exercise 922 (p. 316) 
in Whitworth’s book (loc. cit.). His Q" is the same as Oe ARES 
i.e., aS our at of equation (2). In particular, Whitworth’s equa- 
tions given in his exercises 768 and 769 become, in our notation, 


(eRe am (28) 


and 


y(t) Rear Get ay a". CY) 


The second of these follows easily from the first since () = 


(i > 1) = ( 7) . To prove (28), we use equation (2). This 
mn yp 


gives: 


Be Tether (le-- 
oat Aa a 


*Following Norlund we could also write equation (27) in the form 


n—m 
Se alge) Or Ben: 
m ‘ 

s=i 
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Alternatively, one may use equation (16) and the multinomial 
theorem to prove (28). (See appendix 2). 
Let 

N= 8, 4+ 8, t+--+8, =Gy-Ota,-1 +...4+4,-n, (30) 
with M=a,+G,+..-44,, a,2 0, 
be a partition of n into m positive or zero parts. Then the number 
of distributions of n different things into m different parcels, such 
that s, of the things go into (any) one of the parcels, s, into a 
second parcel, etc., and finally s_ of the objects into a last par- 
cel, is equal to 

n! m! 


n! m ! 


(OU)? (11)%4 3 2 (nly? ’ aha dng hee 


n 


(31) 


In particular, the number of distributions of n different objects 
into m different parcels, with empty parcels allowed, is given by 
the sum of all quantities (31), i.e., m” (Whitworth, Joc. cit., Chap. 
III, proposition XX]). 

Now let a, = 0 (no empty parcels admissible) and consider the 
distribution of n different things into m indifferent parcels. We 
obtain the correct number of arrangements if we divide (31) by m!, 
which gives: 

n! 


(it) fot)" Se tat oe 


(32) 


Hence #” must be the sum of expressions like (32) taken over all 
possible partitions of n into m positive parts, and this is the com- 
binatorial interpretation of equation (12).* 


*Inside each ‘‘box’’ the s; different things may be arranged in sj! 
ways. If these internal differences are taken into account, i.e., if in the 
above statements we replace the word ‘‘parcels’’ by ‘‘groups,’’ the num- 
n! m! 
ber of arrangements becomes —~———_————— instead of (31), and 
Go! a,! ... a,! 

n! ‘ 
Pay ate instead of (32). For example, the number of distributions of 
y! Gag! oe. 

n different things into m indifferent groups, with no empty groups, is 
given by & in!/ay! @,!... a,! , the sum taken over all sets of non-negative 
integers (@1, @o,..., a,) satisfying the two equations (30) with ao = 0. 
, n! [n-1 

The sum is equal to — 
m! \m—-1 


that we have derived Whitworth’s proposition XVII of Chap. II, loc. cit.. 


[cf. equations (34) and (36) below], so 
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The quantity (32) may be written in the form 


n! ' n (ee 
~ = ee ae 
TT (ua {Bat SN a 
p=1 
Now 


n= ¥ (ua,,)3 


yet 


therefore, by a well-known theorem (G. Chrystal, 1952, Volume IJ, 
$12, p. 545) the term in brackets is an integer. On the other hand, 


(u a)! 


Hg! 
(n!) a, 


is equal to = if p = 0 orif o> 0. It is equal to 
a! 


(h 
of aed oy ae 
1 fpv\) _ pv “ 
i Ee) Gea) 
j.e., to an integer, if p > 0 and a,, > 0. Thus (32) is always an 
integer, as of course it should be. 


We shall see that N. Rashevsky’s original problem (to find the 
total number of transformed graphs) consists in determining the 
distribution of indifferent things into different boxes, with certain 
restrictive conditions. We therefore collect some well-known re- 
sults for this type of distributions. 

Given again a partition of the integer n into m parts, like (30), 
then 


oe fae (33) 


is the number of ways of distributing n indifferent objects into m 
different boxes, so that s, objects go into one box, 8, into a sec- 
ond box, etc., and finally s into an mth box. (We may imagine 
putting a label on each box indicating how many objects are to 
occupy it. Since a of the numbers s, are equal to zero, a, equal 
to one, etc., the labels may be permuted in m\/a !a,)..- a,! Ways.) 
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Now let F™'%? be the number of distributions of n indifferent 
objects into m different boxes, with the additional requirement that 
no box shall contain less than qg nor more than q + z — 1 objects. 
Then, using (33), we have 


m!\ 
m3 qyz ss . 34 
re 34) 


! . 
q Gast qt+z-—1 


Here the sum must be taken over all sets of non-negative integers 


ich are solutions of the two equations 
(a,, Cae Bos digewe) which a q 
= : ( = : 
N=a, q+ @,, (0 £1) see ee 34 
Th GG ance a. A Oo a 


By the multinomial theorem (see e.g., Netto, loc. cit., $35, page 
60) the right-hand side of (34), and therefore also ested 8. [- equal 
to the coefficient of #"~%" in the expansion of (1+a2+27+ 
2? +... +a77!)™. We have thus proved Whitworth’s proposition 
XXVIII (doc. ctt., Chap. II). 

On the other hand, let us write ®£'* for the number of composi- 
tions of n into p positive parts whose magnitude shall not exceed 
a. The same number may appear more than once in the composi- 
tion, i.e., these are complexions formed with unlimited repetition 
from the elements 1, 2, 3,..., 4 and with a given sum n (Netto, 
loc. cit., Chap. 5).* Following Netto we shall denote by [®%**] 
the set of all the compositions just mentioned. For example, we 
may write the symbolic equation 


[O27] = {2+ 1h+ {1 +2. 


On the right-hand side the elements of the set are written with 
a ‘+? sign between each of them. Thus ©?2? = 2, 

To find the number of compositions we may first determine the 
partitions of n and then see how many permutations can be ob- 
tained from each partition. Therefore, using again the multinomial 
theorem, we see that ®*** is equal to the coefficient of 2"~? in 
the expansion of (l1+a+a7+...+a%7')* or also to the coef- 
ficient of x” in the expansion of (@ + #7 + #3 +... +4 a%)P.** Itis 


*A more complete notation for ®°~* would be Netto’s cP) ( [n; 
1, 2,..-., &) which, however, is rather cumbersome. : 

**MObius has pointed out the equivalence between certain problems of 
distribution and the composition of number (quoted by Netto, loc. Ctt., 
page 180). See also MacMahon, loc. cit. 
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not difficult to find a more convenient expression for these quanti- 
ties. Let c”” be the coefficient of #” in the expansion of 


eee te. be ett)’. Since 
v-1\w i) 
(lt@+...+2%7')¥ =(1-2")” -—_—— = 
(1 - 2)” 
Fr Ty (- 1) (") eee ghvty 
L=0 v=0 3 i 


rearranging the double sum into a simple summation we find that 
- w\ [uw 1 
De wi (1H +w-1-pv 
Cu ve 2, (") ( w-1 ? Coe 
u=0 


u 

Here o = =| is the largest integer contained in the quotient 
v 

uU . . . . . 

— (o=0 if u < v), i.e., it is determined by 

v 


u=a-v+f with O<B<v-i. 


Using (35) we thus obtain 


pen ive a y coe(f) & & elder ue 1- “A 


=0 (36) 


m-m 
with c= | al. 
2 


and 


(37) 


N. Rashevsky’s equation (30) (1955b) and Whitworth’s propositions 
_ XXV, XXVI, XXVII (Chap. Ill, loc. cit.) are special cases of equa- 
tion (36). 

A direct proof of equation (37), though not a very elegant one, is 
given in Netto’s book (equation (19), Joc. ctt., page 132). A dif- 
ferent proof can be obtained by the EA-method because it is not 
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difficult to show that ‘ - pine a represents the number of com- 
positions of n into p positive parts of unlimited magnitude, the 
compositions being such, however, that specified of the parts 
are > GQ. 

The coefficients c have yet another combinatorial meaning: 
c”™” ig the number of ways of selecting u objects out of w, if 
(a) each object may be taken at most v — 1 times, and (b) the order 
in which the objects are arranged is irrelevant (regular combina- 
tions of order »—1 and of class wu formed from w objects; see 
Netto, loc. cit., $19, p. 31). From this point of view the coef- 
ficients c’*” have been studied in great detail by D. André (1876). 
We shall state some of André’s most important results in terms of 
our numbers 0%, Although Netto discusses regular combinations 
of objects ($$ 19,20) as well as compositions of numbers (Chap. 
5), he seems not to have noticed the connection between them, 
which, however, is only a special case of V. Brun’s general ‘‘duality 
law’’ (Brun, 1927). 

It is clear that ®*°** =0 if 0<n<p or n> pd since in those 
cases there are no compositions with the properties required of 
[®?*“], It is convenient to put ®*=1if p=n=0; Of % =0 if 
p =0,n #0; and also 0°°* = 0 ifn <0,p>0. 

The following theorems are easy to prove from the combinatorial 
meaning of ®/*: 


4, GP = HP” if nen’ =p(a+1). (38) 


We may, in fact, establish a one-to-one correspondence between 
the two sets [oPr“] and [O45 1-,] by associating to each 
composition 


M=A, tAgtecwtA, (1 8A, <a) 


p 


of the first set the composition 
p(X+1)-n=(a+1—-A,)+(A+1—A,)+...+(H41-A,) 


. ya 
which belongs to Cert ee ‘ 


Qa 
o; DP = De BEr ae (39) 


vel 
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Equation (39) is obtained by classifying the compositions according 
to the value of their last part (or any other part in fixed position). 
It is a special case of the more general equation 


ey 


Tr 
OD Pr% x »* Or%. HP-H® (0X7 < p) 


v=rT 


which can be derived in a similar way. 


p 
3. pee). (°) Tastee (40) 


v=0 


For proof see Netto, equation (y), page 1382. (?) pA a aha 
number of compositions of [®**] which contain the part exactly 
v times. 


4. Let y be an integer such that 0 €y<a-1, and let 
a—1)- 
> POLD=A the 


a 
oc aCp— 
Bailye 
oe Ska "0 : (41) 
t=0 p=p—1 


Proof: Using (39) we have 


o o a 
pPs& - p-1,0 
Vy DOs roty * = > Do etasy—y 
t=0 t=0 vel 
It is not hard to see that the double sum can be transformed into a 


simple sum as given by the right-hand side of (41), which is inde- 
pendent of y. 


p(a%—1) oS 
p-1,% 
gt es yi is ge (42) 
b= p—-1 


To prove this, write the expression on the left as a double sum, 


5 3 Dem Bop 


y=0 t=0 


and then use (41). 
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6. The equation 


ap 

Pr% — oP 
ya ag =o (48) 
=p 


follows immediately from (42). See also Netto, loc. cit., page 137. 
%. Combining (41) and (438) we find: 


o 
Pr» = ps 
wy Otaty - 3 


t=0 (44) 
(a-1)- 
Og VS ot: > |e ar 
Oo 
In particular, for 0 =n, p21: 
p vik 
CY as 
Bie ete (45)* 
m= 


8. For fixed p and « the numbers 0/*% increase as n varies from 
p(a + 1) 
ener: [If n increases further the sequence must de- 


crease according to equation (88).] For proof see: André, loc. cit. 
No. 14. Notice that André writes (p, n - p),_, for our ®2°% 

Some more theorems which are less important for our purposes 
can be found in André’s paper together with a ‘‘generalized Pascal 
triangle’’ for the practical calculation of c*” (cf. also Chrystal , 
loc. cit., Vol. I, $15, p. 67). 

We may mention here that the ‘‘inverse equation’’ to (40), namely, 


p 
werent Y* (1)? ( DOr, (46) 


v=0 


is obtained by means of the EA-principle since, for pv 2 1, rae 
is the number of compositions of [®***] which contain the part a in 
v specified positions (the remaining parts being any of the numbers 
1, 2,..., 4). From (46) and (44) one can derive the equation 

p 


yr ogre S é > Un aye = (19-4 


m= | 


which is similar to (45). 


*See also appendix 8. 
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We assume that the reader is familiar with the rules given by 
N. Rashevsky for his transformation (TS? X) (Rashevsky, 1954). 
According to the steps T, and T, of the transformation we have to 
distribute, for fixed values of the parameters n and m, altogether 
n(m —1) subsidiary points, which we may call indifferent objects, 
onto n specialized (distinguishable) points which for clarity we 
shall call different boxes. These latter are in turn subdivided into 
m non-overlapping groups with n, of the boxes belonging to the zth 
group: 


m 
ROMA S15 BL, ms yam (47) 
t=1 


Now, the distribution of our indifferent objects must be such that 
the total number of objects which are placed in the n, boxes of the 
ith group becomes equal to n-7,, empty boxes being allowed. 
Clearly, the number of possible arrangements is not changed if we 
increase the number of objects in each box by one. That is, keep- 
ing the subdivision of the n boxes into m groups as before, distri- 
bute mn indifferent objects into the n boxes, no bow being left 
empty, in such a way that the total number of objects placed in the 
n, boxes of the ith group becomes equal to n. This must be done 
for each possible arrangement of the n boxes into m groups, 1.e., 
for all of the partitions (47), and, thereafter, for each value of m 
from 1 to n. In terms of compositions the procedure may therefore 
be stated as follows. 

From all the compositions of [®**"] select those which have the 
property that their parts can be arranged into m groups, so that 
each part belongs to one and only one group and the sum of all 
parts in each group is equal to n; or, aS we shall say, they can be 
subdivided into m groups of sum n. 

The number of parts in each group and their order within the 
group are not specified. But, of course, the arrangement of the 
parts as a whole is relevant since we are dealing with composi- 
tions. We shall call a composition of [bf>*] admissible if its parts 
can be subdivided into m groups of sum n. Any permutation of an 
admissible composition is also admissible; thus we may equally 
well speak of admissible partitions. For example, ifp=%=n=T7, 
m= 38, 

01=34+24+24+44+7+14+2, OF Theo k, 


is an admissible partition because the part 7 by itself forms a 
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group of sum 7, and the remaining two groups can be taken as 
(8 +4), (2+2+2+4+1); or (4+2+1), (8+2+ 2). However, the 
partition 21=6+3+2+4+2+2+2,0r643 2*, is not admissible 
since it contains only one group of sum 7. 

The number ¥” ™ of admissible compositions is equal to 


! 

p: 
Wm _ ee 3 48 
bas Me {es re Ae =} Sas 


where the sum is taken over all admissible partitions of the number 
mn into p positive parts < a, and the a,’s indicate the multiplicity 
of the parts. 

As an example let us calculate aie If p, denotes the number 
of parts in each group we clearly have only two possibilities: 
(1) p, = 4, p, = 23 (2) p, =p, = 3. In the first case we perform a 
formal multiplication of (41° + 3217 + 21), the aggregate of all 
the partitions of 7 into 4 parts, with (61 + 52 + 43), the aggregate 
of all the partitions of 7 into 2 parts. This gives the following ad- 
missible partitions, all different from each other: 


1) 6414 30 
2) 6321? 120 
8) 62° 1? 60 
4) 54213 120 
5) 53271? 180 
6) 5241 30 
7) 443,35 60 
$).43°9.12 180 
9) 43234 120 


In the second case we have to take the (symbolic) square of 
(517 + 421 + 371 + 32%), i.e., of all the partitions of 7 into 3 parts, 
neglecting numerical coefficients. This gives us ten more admis- 
sible partitions, four of which, however, already appear in our 
previous expansion. Thus we have only six new partitions, namely, 


ny Se 15 
Et} 492.4 90 
12)34 14 15 
13) 3794 15 
14) 53718 60 


15) 3*9*3 60 
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The numbers on the right of each partition indicate how many 
permutations can be formed from it. Adding all of them together we 
find that 


For confirmation of this result we may notice that the total num- 
ber of partitions of 14 into 6 positive parts not exceeding the value 
6 (r<® (f 14; 1, 2,..--, 6) in Netto’s notation] is equal to 16. In 
addition to those we have already found there is, in fact, only one 
more partition, ae is not admissible, namely, 42°. Since this 


gives rise rer) = 6 compositions, we must have 
6,6 _ we,6 _ 
Or, Por =*0 


which is indeed the case because or = 1,161. : 
The number of transformed graphs is given by yr” for fixed m. 
The total number of transformed graphs then becomes 


n 
G = ES gee (49) 


m=1 


this includes the original primordial graph for m = 1. 

We thus have a perfectly well-defined method for calculating G,, 
but it is obvious that the probability of making mistakes increases 
very rapidly even for moderately large values of 1. The computa- 
tions become exceedingly long; they may nevertheless be slightly 
simplified by means of a few theorems which we shall now derive. 

1. First of all we see at once from the definition of ae that 


CE Oe: (50) 


Therefore, using equation (45): 
Gea, (51) 


Unfortunately, however, this upper bound for @, is of no practical 


use whatever, because it turns out that G, increases much more 


slowly with n than n” ax, 


9.lim=t1 


WP = @Pr%, (52) 
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Special cases of this are 


wr als (53) 
sel ae 
#20 = (har) fer. > hs (54) 
3. Ifm=p anda =n 
Pon — Pon _ 
es = ee es (55) 


4, For m = 2 and & =n, we have: 


UP” = OF if p2n+1; 
Wren _ Orn _ peas -—~n, if nis even; (56) 


Stati if nis odd. (57) 


We shall prove equations (56) and (57). The (ordered) partitions of 
[o%”] are classified according to the number of times that they 
contain the part 1. Let 


2h m8, $9, + See St (Rw) ae (58) 


be such a partition, with 


The notation (n-v)x1 means that the part 1 appears exactly 
n—-v times. 

The index v cannot be smaller than 2 because s, <n. Ifv=n 
we obtain the one partition (or composition) 


2n=2+2+4+2+...+2 (n times) 


which is obviously admissible if and only if n is an even number. 
We shall show that any one of the remaining partitions (58), with 
2<v<n-1, is admissible. First of all, s, must be 2 3 if vy <n-1. 
There are now two possibilities: 

a) 8, 2v. Then the part s, together with n — s, (<n-v) of the 
parts 1 gives a group of sum n. 

b) s, <v. Since 


8, +8, +...+ 8 ent+u>d>rv, 
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there must exist an index 7 (2<7<v) with the property that 


Secret toe, but 


Sree Se tet Soe + 8, es 


Here s, is 2 8, and the other s,’s are 7 2; therefore y >3+(r—2)-2 
=2r-1, or 


This shows that r<yv. On the other hand 


See: 
Be tegatana 7S r)-22v. 


L=s 
It follows from this last inequality that 
e=ntv—-pSn, 


Hence, in this case, the parts 8,, 8,,-+--, 8, together with n - « of 
the parts 1 will form a group of sum n (0 €n-eXn-v). 

In a similar manner, although the proof is a little more involved, 
one can show that, forn 2 6, 


Q2n-1 n 
n-l,n—-1 _ n-1i,n—-1 _ ~ a 
ee noo 2. x a € = 4 (n t) (5) (59) 


if n is even 2 6; 
and 


Qn-1 n 
-i1,n—-1 _ @mn-i,n-l _ =. SS a =e 
bh n = OO, n (n Dis Se!) (n - 1) (ie) <4 (60) 
if n is odd? 7. 


5. The equations 


n = —Yn—1 
ELD Ag este (61) 
v=0 
and 
™m 
wert ¥*(-1)? (°) wenn (62) 
yv=O0 


can be derived in the same way as (40) and (46). They hold for 
0 <m &p if we assign the values 


WP = TOF meat, p> leris 1; nl; 


oh Pag eee > aa 
dng 9 for p=em=0, a2l1, 
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Equation (61) is quite useful, especially in conjunction with the 
following theorem. 


6. sf eh - 0, 
if there exists no partition of p into m positive integers all of 
which-are 2 = and <n. 

For if 
i (63) 


is a composition of [W*"] and p; denotes the number of elements 
in the jth group of sum n, we have 


Lepore a tien 2. 


p=. p, tupg + 2.a > pL (64) 


Now let p,, pg, +++» Hp, be the elements of (63) which belong to 
the jth group of sum n: : 


By t Hg tore thy =N- 


Since the p’s are2 1 and <a, 


An immediate consequence of this theorem is that 
My Gaon’ if p< Sm. (65) 


Using this result equation (61) can be written in the form 


K 
mS p sage tag 
Lata = = f- 7 la as n-1 (66) 


with 


(67) 
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As an example we may take p=n=8, m=3, 4, and 5. Then, 
using equation (54), we obtain 


8\ (7 8 8 
wan) (@) + (2) mad + (o)e3, 
Aa OA 8 8 
vit-(5) (i) « G) (i)vez + (0) 
ie f5¥47) 4415 8 
ves = (8) (2) « G)uss+ lyse. 


Here warn is known from equation (59). 

7. The calculation of hea In this case each group of sum 
n must consist of exactly two parts: (n — 1,1), (n - 2,2), etc. Take 
all the compositions of ee increase each part by 1, 
and then insert v times the part 1 and v times the part (n - 1). 
This can be done in 


( hes 4 2 2p\ (av 
V V 2v Vv 
ways (Netto, loc. cit., $23, page 39). Taking v = 0, 1, 2,..-, 4, 


we obtain all the compositions of [V7 Hn-1], Thus we have the 
recurrence equation 


u 
- 2p 2v 2 —v),n—-8 
wpnect =” (Ge) (7) east. (69) 


y=0 


re 
oom 
ow 

+ 


2 
. PM iste; Ae 2 ph, 2) ‘ad 
Since fe =1 for all values of » 70, and ¥7', (1) (ef. 


equation (76) below], again for all values of p, this recurrence 
equation provides a much simpler way to calculate the quantities 
ee than the original equation (48). Moreover, equation (48) 
can be simplified as follows. 


n 
a) If n is an even number and k = as let 


A, tAg tere tA =H with A, 70 (70) 


be a composition of p into k positive or zero parts. 
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Taking 
the parts 1 and (n-1) A, _ times, 
the parts 2 and (n-2) dA, times, 


! n n : 

the parts 6 - and G + ) A, times, 
and finally 

the part - 2A, times, 
we obtain an admissible partition of ke n—1]. Therefore 

(2)! } 
wes he { Sa talbioo igtd) 
Hyn pe Obeikqbarieten ote aw 


(At. FAKED 
A, 20 


Arranging the compositions (70) according to the value of their last 
part and writing the right-hand side of (71) as a sum over partitions, 


a,:O+a,-l+...+@ -p=n, 
. s } (72) 
An +a, +...+a,=k, a,70, 
we find: 
Py he 
fe : 
sg (2 2)! (x - 1)! 2 a, 
2a 2a err ‘ 7 es 
Caastayrcgegy (Ol) “BN Mul) Ee agbee hace! ao (*" 
Vv 
(73) 


For example, assume » = 2. Then, if x 2 2 we have two possibilities: 
Gyeex a1 OTe OP taal and 
@,=k-2, a@,22, a, =0, 


so that equation (73) becomes 


Cn 4! af 4!(x — 1)! 2 
Nan = @pa fx-o+ g]+ 29 lm -2)+ 5| 


Lee = 38(n-2)(n-1)+1  (n even). (74) 


or 
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Obviously this equation holds also for x = 1, i.e., n= 2. In the 
same way we find for all even values of n ? 2: 


if p=3 
went = 5(n - 2)[3n(n - 4) + 14) + 1; (75) 
if p=4 
ws.n-1 = 
4,n 


T(n — 2)[15(n — 8)(n — 6)(n — 4) + 150(n — 6)(n — 4) +295(n — 4) + 79] +1; 
and if p=5 
Y Sele = 3(n — 2)([815(n —10)(n — 8)(n — 6)(n — 4) + 
4,725(n —8)(n — 6)(n —4) + 17,825(n — 6)(n - 4) + 
14,910(n —4) +1,492] + 1. 


It can easily be verified that all these expressions satisfy the re- 
currence equation (69). 
b) If n is odd 


2pjn—1 
+ ee 


(2 1)! x! 
C8 
Yee Le hoe oe I “ 


(dg: Gy9-++s Fy) 


the sum being taken over all partitions (72), where now k = is : : 
In this way we obtain: 
for n= 2 

eX. =3(n-2)(n-1) (n odd); (77) 
for p=3 

i = 5(n -1)[8(n-1)(n- 4) + 8}; 
and for p = 4 
a 


35 (n <1) [3 (n 7) (n — B)(n — 8) + 18(n ~ 5)(n - 8) +17 (n — 3) + AI. 
8. Equation (66) for m = p — 1 gives 
wen i =(n-1) (5). | (78) 


-i,n 
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9. For m = p — 2 we find in the same way, using (74) or (77), 


Oo eee (5) (" r + () [3 (mn - 2)(n - 1) + 8] (79) 


with 5 = 1if nis even, 5 = Oif n is odd, and (f)-0 if p <k. 
We are thus able to calculate at once the values of G, forn <5. 
The results are 


Gy ie aly 
Galas 
G =. &, 
eel a 
G, = 402. 


If n = 6 we obtain Yi'§ from equation (54), 33 from (59), and W5"¢ 


from (75). Inserting all this into (66) we can calculate ae 
Finally, it is found that 


G, = 4,829. 


The results for n = 7,8,9 are shown in the tables below; for com- 
parison we also give the corresponding values of ®"" . 


TABLE 1. n=T7 


m Wnt Om 
1 1 1 
2 1,708 1,709 
3 23,835 26,769 
4 24,430 60,691 
5 3,675 26,769 
6 126 1,709 
1 1 
Total G_ = 53,776 117,649 = 7° 


TABLE 2. n= 8 


G6ebheG6GV3373uNKeewwwwowoaooaoa@s@sS$<S09aS SSS9amaaSaS 


2 6,427 6,427 
3 188,370 193,705 
4 479,571 848,443 
5 146,958 848,443 
6 10,066 193,705 
7 196 6,427 
8 1 1 


Total @, = 831,590 2,097,152 = 87 
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TABLE 3. n=9 


. Tia Om 
1 1 1 
2 24,300 24,301 
" 1,325,917 1,343,521 
4 7,320,168 10,350,421 
5 4,744,740 19,610,233 
6 622,986 10,350,421 
7 23,520 1,343,521 
8 288 24,301 
9 1 1 
Total G, = 14,061,921 43,046,721 = 9° 


An approximate value of G,,. was obtained as follows. Writing a 
complete scheme like (68) we estimate the quantities gt Siar 
along each column by means of the equation : 


weetK+tin—l ~ W2Ht+K,n—-1 2u4+K+i1,n—1 _ 2u+K,n—1 
Lyn yn * iis Din 1. 


We start with x =0 for which the exact value ee is known. 
This procedure gives 


Gio ~ 368 x 10°. (80) 


The values obtained in the same way for n= 8 and n =9 are: 
G, = 851,000, G, ~ 16.9 x 10°. Therefore the right-hand side of 
(80) is probably too high. However, we can certainly say that the 
order of magnitude of G,, is 3 x 10°. 

Equations (49) and (66) give, after a few simple transformations 


ae, od : [5] 
aes n A,n-1l < n r,n-1 
f,.% 2B (7) os Snee es (") De Danas 
A=0 B=0 r=0 L=0 
We may therefore expect that if n is sufficiently large 
nea! 
it nae 5 


Even this, however, is not a very useful upper bound, as can be 
seen from our tables. 

Thus it seems very difficult to survey the behavior of G, as n 
increases. Nevertheless, the few cases we have investigated are 
sufficient to cover the range of values that were of most interest 
to N. Rashevsky in his work on the total number of transformed 
graphs. 
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Appendix 1. We have 


m" noi _(m—-A)" 
2 ae Dee ee 


From the first of Netto’s equations (17) (Joc. cit., page 249) it 
follows that 


n —-1 m-—-1 
m! a A!(m — A)! 


If we insert this into the preceding equation and notice that 


( m ee ; y n-m ( m ): 
m—aA a (m — A) pee wy fp 
v=0O 


we find 


kas Vv 


em “i As DM ry! 


Here the second sum is equal to pet en é 


1 
Appendiz 2. From equations (1) and (16) we obtain 


If 


x Ge (2) 


is a composition of n into u positive parts, we may insert (m — 1) 
times the part zero to obtain a composition into m non-negative 
parts. This can be done in 


ra) aCe) 
m — pL I 
ways (Netto, loc. cit., $23, page 39) for each composition (2) and 


for » = 1,2,...,m. Therefore, the right-hand side of (1) is equal 
to the multinomial expansion of (1+ 1+... +41)", i.e., of m” 
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Appendie 3. Expressing ©" by means of equation (37) we see 
that o <m-1(p 2 1). We may, therefore, write 


apne o\ [(m - p)n-1 
PPm — =1)" SAG tide 
Sapien Sula ease 


The right-hand side can be transformed to 


ae scenall be By see 
2. a A (es | eh) side (a) 8 (). 


Here the inner sum has the value (- 1h aaa (? ie! ‘ 
p _ 


equation (14), page 248). Thus, putting p - A = 8, we obtain 


p p,n es s pn-1—sn (Dae 
Bone Fem OCs) 


m=1 


(Netto, loc. ctt., 


This is equal to n®~* if n21 (Netto, loc. cit., equation (45), 
page 256). 
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